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Set Books 


G. D. Smith. Numerical Solution of Partial Differential Equations (Oxford. 1971). 
H. F. Weinberger, 4 First Course in Partial Differential Equations (Blaisdell. 1965). 


It is essential to have these books; the course is based on them and will not make sense 
without them. They are referred to in the text as 5 and W respectively. 


Unit 1 is based оп W: Chapter I, Sections 1, 2 and 4. 


Conventions 
Before working through this text, make sure you have read 4 Guide to the Course: 
Partial Differential Equations of Applied Mathematics. 
References to Open University courses in mathematics take the form: 
Unit M100 13, Integration II for the Mathematics Foundation Course, 


Unit M201 23, The Wave Equation for the Linear Mathematics Course. 
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1.0 INTRODUCTION 


The idea of wave motion, a periodic phenomenon in time and space. must be familiar 
lo everyone. In many instances it can be observed by the unaided eye as in the case of 
waves on the surface of water or the movement ofa disturbance along a rope after the 
end has been moved suddenly. In other cases, the motion of individual waves is 
impossible to observe directly, as in the case of vibrations of a very taut violin string 
or sound waves in the air. Only after detailed experiment and analysis can it be shown 
that the phenomenon is a type of wave motion. One feature common to all these types 
of wave motion is that energy is transferred through the medium without perma- 
nently displacing it. Even if no medium is present, such as in the case of electromag- 
netic waves in vacuo, energy is still transferred. 


In the physical world as we understand it there are three spatial dimensions and one 
time dimension. Variables describing physical quantities which vary continuously in 
space and time arc, in general, functions of all four of the coordinates necessary to 
specify a point in the four-dimensional vector space making up the space-time 
continuum. 


The use of all four coordinates to analyse the wave motion tends to obscure the basic 
relation between the time and space variations which is characteristic of the motion. 
We therefore adopt the usual approach and investigate firstly one-dimensional wave 
motion in which only one spatial coordinate is used. 


You have met the one-dimensional homogeneous wave equation in Unit M201 23, 
The Wave Equation. To refresh your memory, it is the equation 


where u is a function of two variables and c is a constant. Any function u satisfying 
this equation is a solution of the wave equation with parameter c. 


The question of suitably defining a domain for u and the need for imposing further 
conditions on u before a unique solution can be obtained will be discussed in this unit. 


In Section 1.1, we start by deriving the wave equation for three different physical 
situations in order to illustrate the importance of the equation in applied mathematics 
and to indicate some of the methods which could be used to derive the wave equation 
in other situations not discussed in this text. In Section 1.2 we discuss a method of 
solving the wave equation. 


PDE I.L0-1.1.1 
i.l PHYSICAL EXAMPLES OF THE WAVE EQUATION 


1.1.0 Introduction 


You will recall from Unit M201 23 that the wave equation arises when a mathematical 
modelling process is applied to: 

(a) the longitudinal vibrations of a heavy spring: 

(b) the current and voltage oscillations in an electrical transmission linc. 

There is no need for you to re-read the relevant sections of Unit M201 23 in detail 
but a brief glance at them may refresh your memory and help with what follows in this 
unit. 

In the next three sections we are going to apply the processes of mathematical model- 


ling to a further three physical situations and demonstrate that the wave equation 
provides an approximate mathematical model of the situation in each case. These are: 


(a) the transverse vibrations of a stretched string: 

(b) the longitudinal vibrations of an elastic bar: 

(c) one-dimensional sound waves in a gas. 

The wave equation also arises from other physical situations, many of which require a 
wider and deeper understanding of physics than we assume for this course. We hope 


that the few examples we do provide are sufficient to highlight the importance of the 
wave equation in applied mathematics. 


1.1.1 The Transverse Vibrations of a Stretched String 


Ifa string is stretched and the ends attached to fixed supports A and B, it is intuitively 
obvious that a possible state of the string is that of equilibrium, i.e. each particle of 
the string lies on the straight line between the supports and is in the same position 
at all times. 


To study the non-equilibrium positions, we use a rectangular coordinate system (x, y) 
for the plane so that support A lies at the origin and support B on the x-axis at x = l. 


ya 
j4- 1 = 
——————————— —————— 
A 
B x 


String equilibrium position 


A particle of the string can be displaced from its equilibrium position longitudinal] 

(parallel to the x-axis) and/or transversely (parallel to the y-axis). We consider i b 
transverse displacements. If a segment of string is displaced transversely froi n ita 
equilibrium position it exerts forces on the neighbouring segments of а Sith E: 
side. As a result, these are no longer in equilibrium. They too will move ус 
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and so the effect of a transverse displacement is propagated in both directions along 
the string. 


To look at the situation in more detail we set up a mathematical model of the string. 


Let the string have uniform line density p (i.e. mass p per unit length). Let the tension 
in the string be T and let the transverse displacement of the particle of the string with 
coordinate x be y(x, t) at time г. For small displacements у(х. г) we may assume that T 
is constant. We consider the small element of string between points P and Q with 
coordinates (x, y) and (x + Ax. y + Ay) respectively* at time t. Let the length of this 
element be As and let the tangent to the string make an angle y with the x-axis at the 
point P and ij + Ay at the point Q. We see that 


W(x + Ах.) = rit) + Av(x. t). 


Now we require Newton's Second Law of Motion which you met in Unit M100 31, 
ifferential Equations 11. Briefly the law states that 


FORCE — MASS X ACCELERATION. 


We shall apply Newton's Second Law of Motion to the transverse motion of the 
element of string PQ. 


A FORCE acting in the plane may be represented by a geometric vector F (see Unit 
M100 22, Linear Algebra 1). If i and j are unit vectors in the x- and y-directions 
respectively, then {i,j} is a basis for our space of geometric vectors. We then have 


F = Fi Ej 


where F, and F, are uniquely determined real numbers which we call the x-component 
and the y-component of the force F, respectively. It should be clear that if F represents 
a force whose magnitude is F and whose direction,makes an angle 0 with the positive 
direction of the x-axis when measured counterclockwise then 


F, = F cos 0 
and 
F, = F sin 0. à 


So the y-component of the total force acting on the element of string PQ is 
Тїп + Ay) — Tsin y. 

The Mass of PQ is pAs, and since 
As = length of the chord PQ = Ax sec y 

this mass approximately equals pAx sec y. 


* [n this course, as in the Linear Mathematics Course, we shall use the simplifying convention of denoting 
the image under a function by the symbol Гог the function in circumstances where there is no confusion 
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Finally. the ACCELERATION of PQ in the y-direction at time z is given by the second 
derivative of the function 


t— у(х.) (2 0. 


which is just 


We shall now proceed to the limit as Ax and Ay tend to zero. We have 


, siny Ay)-siny d. i 
lim ——— Aj = di (sin y) 
= cos iy. 

and - 

ü Aplix) _ lim W(x + Ах.) — Wx. 0 

m Ax o Ax 
= a (x. t) 
ёх 
This gives us 
р 
T cos^i ACA (1) 


We know that at any instant of time 


(2) 
Finally, 


соңу = эз лл. = 1 . 
sec" — (1o лапу) [1 + (reo TF 


If the displacements y(x, г) are small compared with the length of the string and y isa 
fairly smooth function of x, then (Cy; Cx)? « 1 and costy = 1. (The symbol « means 
"very much smaller than".) This assumption is equivalent to the assumption that y 
remains small. Equation (2) becomes. approximately, 


or 


where c 
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Equation (3) is the wave equation in one dimension. It describes the transverse motion 
ofa stretched string quite well provided that the displacement у(х. г) is not too large. 
Note that it is a linear differential equation. Thus. if y, and y, are any two solutions, 
then A; y, + 4,5» is a solution for all choices of the constants A, and Ap. 


In addition to satisfying the wave equation, the transverse displacement у(х. г) must 
satisfy various initial and boundary conditions. At the fixed supports A and B on the 
x-axis, where x = 0 and x = l respectively. we have y(0, t) = y(i. () = 0 throughout 
the motion. These conditions are a particular form of boundary conditions, i.e. con- 
ditions at a spatial boundary for all times. In addition, at the time t = 0 we can start 
the motion of the string with a given shape 


y(x.0) = f(x) 


and a given velocity distribution 


dy eo cau 
2100) = g(x). 


These conditions are a particular form of initial conditions. 


We shall see in Section 1.2 that these initial conditions and boundary conditions 
constitute a set of conditions that determines uniquely the solution to the wave 
equation for the string for all future time along the whole length of string. This prob- 
em, formulated mathematically as follows, is the vibrating string problem 


(0<x<l,1>0) 


(<x<l) 


0) 
yh) 20 (r 2 0) 


Find a function y satisfying the equations above and continuous in the domain 
(0, I] x [0, co). 


In W: pages 1 to 5 the vibrating string problem is derived in a more rigorous way and 
in more general circumstances. We feel that Weinberger's treatment is too complex 
and long for inclusion in this course and even then it does not tell the full story. 


SAQ 1 


If y(x) is the angle between the tangent to the curve у(х. to) at some fixed t, and the 
x-axis, and s(x) is the length of the curve from some fixed point to the point (x, у), we 
can define the curvature K(x) of the curve at the point (x, y) by the equation 

vo). 

s'(x) 


Find К in terms of y and у”, where y’ is the derived function of x — у(х. to). 


K(x) = 


(Solution on p. 32.) 
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1.1.2 The Longitudinal Vibrations of an Elastic Bar 


READ W: page 6, line 20 to page 7. line 10. 


Notes 


W: page 6. lines — 13 to —7 | 
We show how Си, x(x. 0) is to be interpreted as the elongation per unit length 
at the point x at time г. 


Consider the element of the rod which lies between the planes x and x + h in 
the equilibrium position, where h is small, as shown in the diagram. 

Equilibrium 
position 


» 


Displaced 
position 


4———» Б таи e 
"m u(x + Л, 0) 


The length of the element in the equilibrium position is 
lop =(xth)—x=h. 

The length of the same element in the displaced position at time t is 
1, 


(u(x + һи) + (x + Л) — (ибх, t) + x] 


u(x + ht) — u(x,t) + h. 
Therefore the elongation of the element per unit length at time г is 


= 1 u(x * h.t) — u(x, 1). 


h h 
We now take the limit as л becomes small to obtain the elongation per unit length 
at point x at time t: 


u(x +h, t) — u(x, t) " Qu 


= (xt) 
ôx 


lim 


h~o h 


We do not propose to justify the statement in W that the tension T depends 
linearly on би/ёх (and not on higher derivatives of u). We shall adopt the view that 


Tix.) = E (x. 1) 
ёх 
is an empirical result (i.c. based on experimental evidence) which is valid provided 
Cu/Cx remains small. 
W: page 6, line —6 


A body force is a force which acts on each element 


| [ of a body from outside the 
body. The usual example is a gravitational force. 
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(ii) W: page 6, line —4 
We can write this equation as 


ди ĉu » хз eu 
Ex t)— Ex Кы!) + T p(x) F(x, dx = [. р(х) az (х, t)dx. 


1 
1 
н 
! р(х)Е(х.)Ах 


D amens 
и e" 


The left-hand side corresponds to the total force (per unit cross-sectional area) 
in the x-direction acting on the portion of the rod between the planes:x = x, 
and x = хз. The right-hand side is the product of mass (per unit cross-sectional 
area) and acceleration integrated over the 'same portion of the rod. It is the 
analogue, for a continuous medium, of 

аи; 


m; — (t) 

D бш? 

for a system of discrete masses m; and corresponding displacements u(t). The 
continuous variable x replaces the discrete variable i as a labelling parameter. 


The equation represents the statement of Newton's Second Law of Motion for 
our particular situation. 


(iv) W: page 6, line —2 
This equation comes from the previous one by means ofthe Fundamental Theorem 
of Calculus (Unit M100 13, Integration II), applied to the function which maps 


x2 o? 
X3 — | das - ғ} 
x et 


We can write the result as 


Note that c is a function which is constant only if the mass density p(x) does not 
vary with x. If F is not the zero function, the above equation is nonhomogeneous: if 


Е:х э 0 O<x<i, 
the equation is homogeneous. 


(v) W: page 7, line I 
This sentence refers to boundary conditions, which were discussed briefly in 
Section 1.1.1. For our problem they are dealt with in SAQ 3 and SAQ 4. 


(vi) W: page 7, line 2 
The equilibrium solution referred to here has a constant tension Tọ associated 
with it. It is a different equilibrium solution from that used in the definition of 
the displacement u(x, г). The latter is the equilibrium solution with zero tension. 


SAQ 2 
Suppose that, in the analysis of the motion of the elastic bar, we replace the relation 
ĉu 
T=E— 
ex 
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where K is constant throughout the motion. Find the equation which replaces the 
nonhomogeneous wave equation 


Is it linear? 


(Solution on p. 32.) 


SAQ 3 


An elastic bar, bounded by the planes x = 0 and x = |, has its ends fixed absolutely 
rigidly to a table top made of incompressible material. What are the boundary con- 
ditions appropriate to the longitudinal vibration problem of the bar? 


(Solution on p. 32.) 


SAQ 4 


An elastic bar. bounded by the planes x = 0 and x = / in the equilibrium position, 
has free ends What are the boundary conditions appropriate to its longitudinal 
vibrations? 


(Solution on p. 33.) 


1.1.3 One-dimensional Sound Waves in a Gas 


The final example of wave motion that we consider is sound waves. Sound waves 
differ from waves along a string in that they are longitudinal waves, as in the vibration 
of the bar studied in the previous section. The molecules of the air move in the direction 
of propagation of the wave, so that there are alternate compressions and rarefactions 
of the gas, and the restoring force responsible for propagating the wave is simply the 
opposition of the gas to compression. 


In this example we shall consider only plane waves, or waves having the same direction 
of propagation everywhere in space; waves inside tubes of uniform cross section will 
usually be plane waves, also sound waves a long way from the source will approximate 
to plane waves. We let x be the coordinate parallel to the direction of propagation. 


To obtain the relevant equations of motion of the gas we first find a more general 
equation of the motion of an incompressible fluid in one dimension and then make 
the approximation that the motion departs only slightly from the equilibrium motion. 


In order to describe the one-dimensional motion of a small disturbance (sound wave) 
in a gas we shall require, in addition to Newton's Second Law of Motion. an important 
physical law which we have not yet encountered. This is the Law of Conservation of 
Mass, which states that 


the total mass of a quantity of material remains constant in time, 
however it becomes distributed in space. 


Let x = x, be a fixed plane and (for x; > x,) denote the total mass of gas between 
the planes x = x, and x = x, at time t by M(x,, х. t), as shown in the figure. 


Let u(x, t) be the velocity in the x-direction at time г of the gas passing through the 
plane x = constant. 


If p(x, 0) is the density (mass per unit volume), then gas enters the cylinder bounded by 
the planes x = x, and x = x; (x2 > xj) at the rate Au(x,, р(х. t) and leaves it at 
the rate Au(x5. r)p(x5, t), where A is the cross-sectional area of the cylinder. The net 
rate of increase of mass in the cylinder is 6M/@t(x,, x2, t). Therefore 

a 


M 
oo 1x, 1) = Alu(x,, р(х, 0) — u(x, plx, 0] Vx > x, 
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pix,. 0 (хз. 1) 


P(X2, t) 


u(x, 0) 


From the definitions of M and p we have 
A- (x, X. 1) = Ap(x. 0). 
N 


Now we differentiate the expressions for ёМ/бг and 2М/ёх partially with respect to x 
and ¢ respectively, thus obtaining 


OM [up] 
= = - А 2 
Cxét ёх 
апа 
ё 
= 458. 
ct 


Assuming that M is sufficiently well-behaved*, we obtain 
PM М 
ET 


CXCU 


which leads to 


ep 
et 
or 
cu ep ёр 
pat eS: 
ex ex ё 


This is the mass-conservation equation, which represents mathematically the appli- 
cation to our problem of the Law of Conservation of Mass. 


We now consider the application of Newton’s Second Law of Motion to the gas in 
the cylinder x, < x < x». 


To obtain the correct expression for the acceleration we introduce the concept of 
differentiation following the motion. Let f(x, t) be some property of the gas, such as 
velocity or density. Then éf/¢t represents the rate of change at a fixed plane along 
the x-axis. The rate of change of f for a given element of the gas, whose position x 
varies with t, is given by the total derivative df/dt. We call this rate of change the 
derivative following the motion. 


For a general variation of x and г we have the first-order Taylor approximation 


inp лх+ м 
ex et 


(Unit M201 14, Bilinear and Quadratic Forms). Thus 


* See Unit M201 14, Bilinear and Quadratic Forms, Section 14.2.3. 
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ef 


Af Ê — + 
ex м б 


rom 


and in the limit as Af. Ax and At approach zero 


df Ordx Lf 
d ct 


where u gives the velocity. 


This result is a particular example of the chain rule for composite functions of several 
variables*, which is a generalization of the chain rule which you met in Unit M100 12, 
Differentiation I. 

Suppose 


ух, 


апа * 

ep. 4) — x. 

1:0р,9) Y. 

Cp. g) — 2. 
where (р, 4) belongs to the domains of č, y, C. Then the chain rule tells us that 
Qu ez 


z Op 


and 
Quz. 
д: да 


Here би/др means the derivative at p = p of the function 


pr м(ё&(р, qun. a.p.) (BE P); 
ди/дх means the derivative at X = x of the function 
х w(X, y, 2) (Xe X); 
and so on. (P and X are the relevant domains.) 


The rule can be generalized to any number of domains and any number of com- 
positions of functions. So, if 


0:11 c— x, 
Buy 
ytez 
and 
SAX, yz) e vs 
then we have 


dv _ Ov dx Qv dy | Ov dz 
dt Ox dt. дуй — Oz dt 


judiciously choosing d or ĉ according to whether the relevant function has a single 
or multi-variable domain. 
The acceleration of the element of gas which, at time t, occupies the disc (of the 


cylinder) at x is given by differentiating the velocity following the motion. Thus we 
obtain 


* For a derivation of the chain rule for a function of several variables see a i 
text о T 
W. Kaplan, Advanced Calculus (Addison-Wesley, 1962). кинин 
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a(x.t) 2 u(x. pec (х.ї) + с (х, 0). 
Cx ct 


To apply Newton’s Second Law of Motion we must specify the forces on the region of 
gas. The total force in the x-direction at time r on the cylinder of gas bounded by the 
planes x = x, and x = xj (x, < xj) is composed of two different Lypes of force. The 
first is due to the pressures (i.e. forces per unit area) across the ends of the cylinder. 
both of which are directed inwards: the second is the body force represented by the 
integral of the product of the density р(х. г) and the body force per unit mass Р(х, 1) 
over the whole of the cylinder. So the total force is given by 


ху 
Ap(xi.t) — Ap(xs. t) + af р(х. OF (x, t)dx, 
x 
where p(x, t) is the pressure across the plane х = constant at time t. Newton's Second 
Law of Motion, applied to the total system. gives 


ax. 


pix. t) — р(х, t) + | pix DF, Dax 


х2 д, д, 
= || p(x, of ts t) C(x, ) + Six, ole 
A Ox e 


after substituting the expression for the acceleration obtained previously, and dividing 
throughout by A. Partial differentiation of this equation with respect to x; gives 


Aen. + pixa, OF(xs t) = С + Fenn} 


or, in terms of functions, 
Qu фр 
= + puli = pF. 
Pa 4 тах? 
The two equations which we have derived here will be used іп the next reading passage 
to obtain the equation for sound waves in a gas. 


READ W: page 7, line 11 to page 8, line 7. 


Notes 


(i) W: page 7, lines — 17 and — 16 

The pressure at a point in a gas depends on the density and the temperature. In 
general the temperature may vary independently of the density. If the gas is 
vibrating so rapidly that there is little tinfe for much heat to be transmitted 
through the medium in one direction before the flow of heat is reversed, then the 
system is adiabatic. In these circumstances and with an initially uniform tempera- 
ture distribution it can be shown that the pressure is completely determined by 
the density. (Note how Weinberger expresses this relationship by the applied 
mathematician’s traditional abuse of notation: p = p(p).) 


For an ideal gas the relationship between p and p in the adiabatic situation is 
given by , 


р(х, t) = afp(x.n)*, 


where g is a positive constant, which may be determined from a knowledge of the 
pressure and density at some (х, 1), and y is a constant, greater than 1, which 
depends on the molecular structure of the gas and certain other factors. 


(ii) W: page 7, line —12 
The term p/p is the ratio of two quantities of the same physical type: both are 
densities. It is therefore independent of the size of the unit we choose for the 
measurement of density. Such quantities are dimensionless—they are just 
numbers. Clearly p'(p)/p'(po) is also dimensionless. It is not, however, quite so 
obvious that u/[p'(pg)]* is dimensionless; we shall show that it is. 


(iii) 
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Let M. L and T denote the fundamental physical quantities mass, length and 
time respectively. Velocity has the dimensions LT ! and force has the dimensions 
MLT~2, as we saw in Unit M100 3, Operations and Morphisms. Pressure is a 
force per unit area and so 


dim (force) MLT^? cdm 
күйөр Л ш ыш TT s 
dim ( pressure) dim (area) Le 


where dim is the function which maps a given physical quantity to its dimensions. 
Similarly, 


dim (density) = 


Thus dim maps 
ик Т^! 
and 
HM ML^'TUUMLO = L?T~?; 


so that 
(z) ET: 


So we see that и and [ p(poJ]! have the same dimensions and therefore u/[ por 
is dimensionless. 


It is meaningful in physics to say only that dimensionless quantities are small 
(or large). The sizes of other quantities can be changed arbitrarily by changing 
the units in terms of which we measure them. 

W: page 7, lines —8 and —7 

Equations (1.14) are derived from Equations (1.13) as follows. 


Let 


and assume that A(x, t) and f(x, г) (which are dimensionless) are small compared 
with 1. Substituting for p and u in the first of Equations (1.13), we obtain 


РД 
) 20) — 
Pov P Ро ax 


The second and third terms are small and may be neglected; this leads to 


г ĉi op 
(00) == + ро = 0. 
Pov (Ро) x C POR 


Expressed in terms of p and u, this equation becomes 


ĉu Ср 
= +5 =0, 
бозу et 


which is the first of Equations (1.14). 


The second of Equations (1.13) may be approximated similarly by 


Ou ёр 
nz дешсе eig је 
Po E ix Po 
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By the chain rule for differentiation, we have 


Op _ ёр, 
AC (p) == 


Put 
р(0) = р(ро)(1 + о), 


where х is assumed small compared with 1. Neglecting a, this leads to 
ёи 1. \op 
Poa, + Plo) = PoF. 
This is the second of Equations (1.14). 
(iv) W: page 7, line —4 


We need to assume that 


д2 


9p p. 
Otx — Oxót 


(v) W: page 8, line 1 
We need to assume that 
Cu Pu 


Oxt — tx 


SAQ 5 


In theadiabatic sound wave problem just discussed, find the forms of the wave equation 
obeyed by u and p if the gas is an ideal gas (i.e. p = тр?) and the body force is given by 


F(x,t) = Fo sin ax cos bt. 


(Solution on p. 33.) 


SAQ 6 


In the one-dimensional gas flow problem, if the velocity u(x, г) is independent of the 
time coordinate t, is the acceleration necessarily zero? 


(Solution on p. 34.) 


SAQ 7 


For a steady one-dimensional gas flow (i.e. u(x, t), р(х, t) and p(x, t) are independent 
of time) in which the body force is zero, prove that 


ри = К 
апа 
р = – Ки + L, 
where К and L are constant functions. 


(Solution on p. 34.) 
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1.22 SOLVING THE HOMOGENEOUS WAVE EQUATION 


1.2.0 Introduction 
We now show how to solve the wave equation in terms of the initial conditions over 
the whole future domain for two different types of boundary conditions. 
If 
thea =O (270) 


for a = 0 or а = I. then we have a fixed boundary condition at x = a. If 


xa 


for a = 0 ora = I. then we have a free boundary condition at x = а. We have already 
met physical situations from which these terms arise, in Section 1.1. 


We shall first solve the wave equation 


3 (О <х< [г> 0) 
01? 


subject to the initial conditions 
ulizo = SQ) (0€ x&€I) 


and 


Са 
using d'Alembert's method of solution which was treated briefly in Unit M201 23. 


Note that, to avoid confusion, we are restricting и to be a variable; we shall use v for 
the function 


(x, t) н 


to preserve a distinction between a function and its value. This will be useful when we 
wish to use another set of coordinates to define a function 


мё, n — u. 


1.2. The Fixed Boundary Problem 
READ W: Section 2, page 8 to page 15, line 6. 


The first part of this reading passage (as far as W: page 10, line — 5) is d'Alembert's 
solution of the wave equation with given initial conditions. The remainder of the 
passage determines the solution of the problem posed by including fixed boundary 
conditions. 


Notes 


(i) W: page 9, line 8 to Equation (2.3) 
The change of variables in this passage is most easily followed if we carefully 
distinguish the functions involved. Let 


a(x, i) = & 
and 


B(x.) = л; 


PDE 1.2.1 


we have already. in the Introduction, defined the functions 
v(x. t) —> u 

and 
w(ë, y) u. 

Thus 
wi(alx, t). f(x. 1) — u. 


Applying the chain rule for a function of several variables to the composite 
function w o (x, 3). we obtain 


Suppose now that the functions g and /! are given by 
(x, t) — x + ct, 
Bux )—x-ct 


then 


So 

Ou | Ou. 
Ox дё д 

A further application of the chain rule yields 

u _ O(Qu/Ox) _ 

дх дх 


Similarly, we may obtain 


ё?и Ou Р 
Е с2—— = 0 
e ёх? 
becomes 
2 
aa u - 
ecón 


Since с # 0 (by assumption) we obtain, expressing u in terms of the function w, 


ary ] 
а?у (aec 


been 


Integration with respect to č yields 
éw 
= (ën) = rin), 
ĉn 


where r is an arbitrary function. Further integration, this time with respect to jj, 
yields 


vibus Í наў + ple). 


a 


(i) 


PDE 1.2.1 


where а is a constant and p is another arbitrary function. If q denotes an integral 
of r, we obtain 


u = (Ё, n) 
= рё) + gin) 
where p and q are arbitrary functions. In terms of x and t we have, therefore, 
u = p(x + ct) + q(x — ct). 


W: page 9. line —4 to page 10. line 2 
The travelling wave nature of the functions 


(x, 0) э p(x + ct) 
and 
(x. 1) э g(x — ct) 


is illustrated by the following diagrams. 


^ 
Time t 
"n ad pix + et) 
1 1 
1 1 * 
1 i 
1 Н 
i 1 
1 1 
i 1 
i 1 
1 П 
1 1 
1 i 
1 1 
! i 
i i 
A H i 
Н ] 
1 ‘ 1 
1 Timer + Ar H 
i t 
i 1 
1 Н 
1 1 
i 1 
! „р(х + сч + Ad] 
"i 
i 
1 
1 
1 
1 
1 
1 
1 
1 
1 
S » 
4------ > 4------: > x 
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^ 
Time t 
g(x — ct) 
1 1 » 
1 1 х. 
i i 
1 Н 
i i 
i i 
Н Н 
i i 
i i 
1 1 
! H 
1 H 
Н i 
^ Н Time t + At H 

| 
i Н 
i i 
i Н 
i Н 
i Н 
Н Н 
i i 
Н 1 q[x — c(t + Ad] 
Ц 
1 
Ц 
ү 
I 
f 
E 
1 

— ! > 
а==2=25=ы » acne ducas. > x 

c^t cát 


(ii) W: page 10, lines — 10 to —8 
The inequalities 


OSx+er<l 


< 
O0<x-a<! 


IN 


are, when expressed singly, 

x+ct20, 

XI. d 
0. 
1. 


к= 


AN 


x= 


The third of these is equivalent to 


x 
'©—. 
c 


The second is equivalent to 


We have the basic restriction 
120 


as the wave equation is only being considered in the domain 1 2 0. Since c is 
defined to be positive and г > 0, it follows that 


ct 2 0. 


(iv 


кә 
о 
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Thus 

x + ct 2 0 is a consequence of x — ct > 0, 
and 

x — ci &l is a consequence of x + ct <S l. 


Therefore all four inequalities 0 < x + ct < landO < x — ci < lare satisfied if 


x 
12008 and г 
Т 


The triangular region 


In this region the solution of the wave equation is given by 


х+а 


А 1 
u = хс) + fix — ct) + | g(X)dx 


xe 
and is thus dependent ONLY on the initial conditions (the conditions at t = 0 for u 
and (u/ét for 0 € x < 1) and Nor on the boundary conditions (the conditions on, 
for instance, u at the ends of the string for all t 2 0). Indeed, as far as the solution 
in the triangular region is concerned, the values of v(0, г) and v(l,r) could be 
given by any functions consistent with the initial conditions, and the solution 
in this region would be unchanged. 


We are interested in the solution over the whole domain г > 0,0 < x < I and 
not just in the triangular region of the diagram. In determining the solution over 
this larger domain we shall find that the boundary conditions are needed. 


W: page 11, lines 12 to 16 
We wish to obtain и over the domain defined by the inequalities 


Oxx&l, 120. 
These inequalities are equivalent to 


0<&+л<&<2, 2-2 0. 


PDE 1.2.1 


[UN gest 


(0, 0) 


X 


It is clear from the diagram that this domain is contained in the domain defined 
by the inequalities 


ns i20. 


Thus, a knowledge of p(Z) for ë > 0 and g(y) for y < l is more than sufficient for 
a complete solution to our problem. 


The general results obtainable from Equations (2.9) and (2.10) ( W: page 11) are: 


рё) = —q(2nl — č) (2n — 1) € č € 2nl, 
p(é) = p(é — 2nl) 2nl S Ë € (2n + 1). 
qn = qi + 2nl) —2nl S y S —(2n — 11, 
а) = —p(—2nl — n) —(2n + D)! € y € —2nl, 


for n = 1.2,3,.... 


These equations enable us to find p and q throughout their respective domains 
€ > 0 and y < l, given p(é) and qj) defined only over 0 < č < LO < у < I by 
Equations (2.5) and (2.6) (W: page 10). 


It should be emphasized that the fixed boundary conditions 
ц.о = Uap = 0. (20 


are responsible for giving us these particular results for extending the solution 
into the whole of the required domain. 
] 
SAQ 8 
W: page 17, Exercise 3 
Assume that the ends remain fixed. 


(Solution on p. 34.) 
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SAQ 9 


For the initial-boundary value problem 


2 ot 
E т (exeLht 0). 
CU cx 
„ЖА 
ilo = asin т (<x <I) 
2 = ИХ 
aa = —bsin — (0€ x € I). 
C [i-o l 
MES" MES (t 2 0). 


where a and Р are constants, show that 


[PLI n20,12...., 


т l anei nac + ni 
t, ——|tan^!-— — 
"axo lb с 


Interpret this result physically. 


where 


(Solution on p. 35.) 


So far we have considered examples of the wave equation problem in which initial 
are given. It is possible to have other conditions: for 


T1 cu 
conditions t|,-9 and — 
CU ao 


А ; к eu " ; 3 
illustrative purposes we сһооѕс as given u|,-, and alo i.e. the variables are specified 
á 


x 
on the line г = x instead of the line г = 0. Such conditions could be imposed, in the 
vibrating string problem, by a "position- and velocity-setting device" which travels 
along the x-axis at unit velocity. In the following example the wave velocity c is 
assumed to be smaller than that of the position- and velocity-setting device. i.c. с < 1. 
Example 

W: page 17, Exercise 10 

Solution 


The general solution to the wave equation is: 


и = p(x + ct) + g(x — ct) (0&x&l, г> x) 


Since 
thes (0€ x € 1), 
we have 
P(x. + o) + (0. 0) = ftx) (OS x <1). 
Since 
Z| =0 (0&€ x € 1), 
we have 


єр(х(1 + ©) — eq(x(1 — с) 20 (0€ x € 1), 
using the chain rule. 
Integration of the last equation with respect to x yields 


Ie +o) Cit 
I+ I 


where K is a constant. 


24 
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Solving the two linear equations in p(x(1 + с) and q(x(1 — c)), we obtain 
Kl- 
p(x(t + с) = 41 + oft +50 (0<х< 1), 


апа 


ха - о) = 0 - дл) -ŽE @<х<1). 


We may choose K = 0, since p and q will be added to obtain u and the constant terms 
will cancel. Therefore 


pe) = 31 + ar) O<E<1+0) 


and 


The boundary condition 


И|х=о = 0 (t 2 0) 
yields 
pct) + q(—ct) = 0 (1 20): 
that is, 
90) = —-p(-m (n < 0). 
The boundary condition 
Ulva, = 0 (rz 1) 
yields 
pl + ct) -q(1 с) = 0 п 2 1); 
that is, 


péje-qQ-à (E2140). 
The two equations 

qon = —р(—1}) (1 < 0) 

р) = -42 - č) (2 
are sufficient to extend p(é) and q(y) as far as necessary, because the inequalities 

tx, O<x<l, 


which define the domain over which the solution is required, are equivalent to the 
following inequalities: 


(1(—-с)ё>(1+с), OS E+ <2. 
The region 
n&l-c. €20, 


over which we can determine p(é) and q(1), encloses the former domain and so the 
solution can be found. This is illustrated in the following diagram. 
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0-c61-c0) y=l—e 


(0, 0) 


(Y 


1.2.2 The Free Boundary Problem 


So far we have concentrated our attention on fixed boundary conditions for the wave 
equation. In the next reading passage you will see how the combination of initial 
conditions and fixed boundary conditions can be regarded as a wholly initial value 
problem. 


READ W: Section 4. page 21 to page 22, line —2. 


In our derivation of the vibrating string problem (Section 1.1.1). if the end x = 0 of 
the string is constrained only in the x-direction but is unconstrained in the y-direction 
(for example, if the end of the string is attached to a massless ring which slides on a 
rod in the y-direction), then there is a free boundary condition at that end. Since there 
is no force on the string in the y-direction. the end of the string must remain perpen- 
dicular to the rod: 

ey 

(0,0) = 0. 


ёх 


In general a free boundary condition at х = а is defined by 


SAQ 10 


A string of line density p and tension T; which in equilibrium lies along the x-axis 
with end-points x = 0 and x = I, has a small bead of mass m attached to the end of the 
string at x = 0. The bead is constrained to move along the line x = 0 but its y 
coordinate is not constrained. What is the appropriate boundary condition at x = 0 
for small transverse vibrations of the string? | 


(HINT: Use Newton's Second Law of Motion for the bead.) 
(Solution on p. 36.) 
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SAQ II 


A string of line density p~ lies along the negative x-axis and another string of line 
density p* lies along the positive x-axis. They are joined at the point x = 0. The 
tension in the strings is T. If y7 and y* are the functions representing the transverse 
displacements on the negative and positive x-axes respectively, what are the boundary 
conditions applying at the point x — 0 when the junction is unconstrained in the 
y-direction? 


(Solution on p. 37.) 
In the next reading passage, we shall see how free boundary conditions can also be 


incorporated in the initial conditions in an appropriate way, enabling us to find 
solutions. 


READ W: page 23, line 13 to page 24, line 12 (end of the section), omitting lines —6 to 
—3 on page 23. 


The omitted paragraph relates to characteristics which we shall discuss in Unit 2, 
Classification and Characteristics. 


General Comment 


The type of boundary condition determines the choice of an odd or even extension of 
the functions f and g (the initial conditions). We use an odd extension for a fixed 
boundary condition and an even extension for a free boundary condition. 


Example 
W: page 24, Exercise 8 
Solution 
The function 
x — x1- х)? (хє А) 


is not even about either of the points х = 0 ог x = 1. which is what we require for 
the extensions of f and g, since we have free boundary conditions at both x = 0 and 
xal 


The extension for f given by 
f(x) = (x 2 пр — (x np ((<&х<п+1,п=0,+1,+2,...) 


is even about both x = 0 and x = I. This enables us to extend the domain of f to 
the whole real line. 


^ 


fix) 


-1 0 1 2 3 x 


The function 
x—l (xe К) 
is even about both x = 0 and x = 1 and we shall therefore define g in that way for all 
хє К. We have the general solution (Equation (2.16) in W: page 13) 
xt 


u(x,t) = SEG + 9 + Дх — 0] + if g(X)ds, 


х= 


PDE 1.2.2 


since c = 1. Substituting g(X) = 1 and performing the integration, we obtain 


u(x.t) = 30 f(x +) + f(x — 0 2t. 


ЕТЕ 


Now4 < 2 <5 and -3< -7 < —2; 


and so 


therefore 


and 


13-0107 


using the formula for extending f. This produces 


Їй 


Note that, for fixed x and increasing 1, the terms f(x + t) and f(x — t) are bounded 
in magnitude but the integral of g increases indefinitely. This looks unusual but is 
due to the facts (if we are considering a vibrating string, say) that both ends of the 
string are free, and the whole string is given an initial velocity of unity in the positive 


u-direction 


L (x,0) = gi) = 1. 
et 


The whole string will continue to move with this transverse velocity since neither end 
is constrained transversely. Superimposed on this constant velocity is the wave 


motion due to the initial displacement 
u(x, 0) = f(x). 

SAQ 12 

W: page 24, Exercise 1 

(Solution on p. 38.) 


SAQ 13 

Find u for x e (0, 1], г > 0 such that 
22 32 * 
ae es (0<x<1, t>0), 
дг? дх? 
ul,=o = sininx (0€ x <1), 
ёи 
= = sing €x< 
Blo sin 5лх (0€ x 1) 
ц\,-о = 0 t > 0) 
ди 
Z| =0 > 
Ax lees (t 2 0) 


{Solution on p. 38.) 
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SAQ 14 


In a vibrating string problem for a string of length I. there is a free boundary condition 
at the end x = 0 and a fixed boundary condition at the end x = I. When the initial 
condition functions 


Гоха иро  (xe[0.7]) 


апа 


(x € (0. 1) 


are extended to help solve the problem, what is the period of the extended functions? 


(Solution on p. 38.) 
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1.5 SUMMARY 


The ware equation 


has been shown to describe approximately: 


(a) transverse displacements of a stretched string: 
(b) longitudinal vibrations in an elastic bar; and 


(c) one-dimensional sound waves in a gas. 


The homogeneous wave equation оп [0,/] x (0, zc) has бесп solved over the whole 
future domain (r € [0. < )) in terms of the initial conditions for x e [0. П and fixed or 
free boundary conditions. The general solution is 


xe 


Afix + e) + f(x — e0] + x] g(X)dX, 


n 


where 
fix ulzo (x e [0. 1]). 


and 


px (x e [0, 1), 
Oeo 


and the domains of f and g are extended to the whole real line so that they are odd 
about an end point with a fixed boundary condition and even about an end point 
with a free boundary condition. 
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1.4 FURTHER SELF-ASSESSMENT QUESTIONS 


SAQ 15 


Which of the following statements, about the solution of the homogeneous wave 
equation for a finite string, is/are true? 


A At cach point the solution depends on the initial conditions and the boundary 
conditions. 


B At each point the solution depends only on the boundary conditions. 


C On a certain domain the solution depends only on the initial conditions, but 
elsewhere it depends on both the initial and boundary conditions. 


D On a certain domain the solution depends only on the boundary conditions, 
but elsewhere it depends on both the initial and boundary conditions. : 


E А1 each point the solution depends only on the initial conditions. 


(Solution on p. 39.) 


SAQ 16 


In solving the homogeneous wave equation in the domain 0 € x € 1, ¢ > 0 with 
initial conditions 


g(x) = sin zx ( 
and boundary conditions 
Ш|к=о = 0 ( 2 0), 


find the appropriate extensions of the functions / and g. 


(Solution on p. 39.) 


PDE 1 SAQ 1-3 
1.5 SOLUTIONS TO SELF-ASSESSMENT QUESTIONS 


Solution to SAQ 1 
We know that 
tan y = y’, 


from the definition of the tangent. Differentiation (using the composite function rule) 
leads to 


вес” = y". 


We obtain s' by using the approximate expression 


— = sec (x) 
Proceeding to the limit as Ax becomes small, we obtain 
s = sec y. < 
Therefore 
kt 
5 
OY соз 
7 secy 
= PET teen 
(1 + tan^y)! 
= EE Gm 
ET 
Therefore 
K E 


TOF 


Solution to SAQ 2 


If we apply Newton's Second Law of Motion to that portion of the rod between the 
planes x = x, and x = x; for any x, < xz, we obtain 


, Ou ди la ди „| ди R 
E3020 + SESS - Ex oen - ESSI 


= f° ote] 6.0 — Fenn fas 
А ag ) = Е(х, г) |dx (г > 0). 


Partial differentiation of this equation with respect to x, yields 


Ou ди д?и Qu 
E— + 2K 
ёх? dx ax? Ра? 


pF. 


This is a non-linear equation, owing to the second term on the left-hand side 
Solution to SAQ 3 


Since the table top is incompressible the ends of the bar will always be in the positions 
they occupied with the bar in equilibrium; thus the displacement will always be zero 
at these points. The boundary conditions on the displacement function и are then 


u(0, t) = u(l t) = 0 u 2 0). 


Solution to SAQ 4 


The ends of the bar are free, so that the tension at each end is zero: 


TQ.) = Tl.)-20  ((20. 


Since 


we sec that 


Qu би 
etat was == > 
ax (0, г) x (i,t) 20 (t 2 0) 


are the appropriate boundary conditions. 


Solution to SAQ 5 

We use the equations on W: pages 7 and 8: 
Uu E Ou дЕ 
дг? Ox? д 

апа 
02р 
дг? 


where с? = р(ро). In our problem 


p = ap”. 
Therefore 
dP Z yapt! 
m* yap?" |. 
Therefore 
c? = р(ро) 
ets, 
Po 
where 
Po = %Po 
Also 
F(x, t) = Fosin ax cos bt; 
therefore 
ok = —bFysin ax sin bt 
д 
апа 


дЕ 
“oz = —aPoF cos ax cos bt. 


The required equations therefore are: 


32 2? 
eu 7. ou ^ ^ 
= |260 = = —bFpsin ax sin bt 
po] ex” 


and 


К, 
Be | f£ = —apygFocos ax cos Ы. 


PDE 1$АО 4-5 
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Solution to SAQ 6 


The gas acceleration a(x, t) is given by 
ĉu ĉu 
a(x. t) = u(x t) z- (x, t) + = (x, t). 
ex ct 


If u(x, t) is independent of ¢ then 


so that 


The acceleration is independent of time but is zero only if u is a constant function. 


An element of gas moves through space and can change its velocity in doing so even 
though the velocity at each point in space may stay constant. For example, consider 
gas flow in a cylinder in which the density varies along the axis but does not change 
with time. 

Solution to SAQ 7 


Using Equations (1.13) on W: page 7 with F = 0, putting time derivatives equal to 
zero and using a prime (') for 0/dx we obtain 


ри + up = 0 

апа 
puw + р = 0. 

The first of these equations can be written as 
(pu) = 0, 

which gives 
pu — K, 

where K is a constant. 

Substituting this result into the second equation, we obtain 
Kw +p = 0. 

Integration yields 
Ku+p=L, 

where L is a constant, i.e. 


p= —Ku+L. 


Solution to SAQ 8 


We have 
OU д?и 0 " 
КТЫ rm (0cx«l r»0), 


ullo = f(x)-sinzx (0<x <1), 


eu 

= = g(x) = <x< 
EM g(x) =0 (0<x <1), 
ці, о = „1 = 0 (t 2 0). 
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We can extend f and g as odd functions about x = 0 and x = 1 with domain the whole 
real line by defining 


Гбх) = ѕіплх — (xeR) 


апа 
g(x) = 0 (хє К). 
This satisfies the boundary conditions. Using the general solution 
] pete 
u = A[f(x + ct) + f(x — en] + zÍ g(X)dx 


we obtain 

u = Hsin x(x + ct) + sin n(x — сї)] 

= sin zx cos лс. 

(We have used the trigonometric identity 

sin(4 + B) + sin(A — B) = 2sin A cos В.) 
The solution to this problem is an example of a standing wave. The outline of the 
string at any time is a half-period sine wave and all particles of the string oscillate in 
phase with the frequency 5c, but with different amplitudes. Roughly speaking, it is 
the two-dimensional side view you would get when watching a skipping rope being 
swung round in circles between two children, each holding an end, moving their hands 
just sufficiently to keep the rope in motion. 
Solution to SAQ 9 
We use the natural extensions for f and g, which are odd about 0 and 1, for the fixed 
boundary conditions: 

О >; 
f(x) = asin T (xe А), 


and 


g(x) = -b sin (x e R). 


The general solution 


u = fx + ct) f(x —c)) + xf g(x)dx 


produces 
al. л(х+с) . n(x—ct) bopee cxx - 
u = zf sin +sin sin — dX 
2 1 1 2с) 1 
.ORX nct Ib . mx . net 
asin — cos sin — sin 
1 пс I 1 
. 7X nct 16. nct 
= sin a cos — —sin x 
I { пс 1 


(We have used the trigonometric identities 
sin(A + B) + sin(A — B) = 2sin А cos B, 
cos(A — B) — cos(A + B) = 2sin Asin B.) 


Whenever the time dependent factor is zero, that is, for values of t satisfying 


и is zero for all x e [0,7]. Therefore 


PDE I SAQ 9-10 


ul, = 9 xe[0.I] 
at time г, given by 
1 лас nl 
t = лап" |+ (20.142... 
пс ib c 
The time dependent factor in u can be expressed (after a considerable amount of 
manipulation) as 


nct db . nct P lb\? . [net 
acos—- = —sin — = — Ја? + |—| sin|— — ф 
1 ле l [rs 1 


лас > . 
where tan ф = E Thus u has the form of a standing wave (see solution to SAQ 8) 


which occupies the equilibrium position at times t = t, (n = 0, 1. 2,...). 


Solution to SAQ 10 


Let у(х, t) be the displacement at time t of the element of string at x. The acceleration 
of the bead is Йй 


— (0,0), 


Y 


А 


t» 


t 


its mass is m and the y-component of the force acting on it is Tsin (г) where y(t) is 
the angle at time г between the tangent to the string at x = 0 and the x-axis. Newton's 
Second Law of Motion gives 


Tsin p(t) = mz. D). 


Now 


sin y 


cos y {ап 


tan ij 
sec y 


tan y 


Jl + tan?y 


and 


Thus 


A2 


ey 
=m J 


Ay H a 
r«| 2o] Е 
сх 


(0, t). 


we obtain, іп this approximation, the linear boundary condition 


2 


[^ 


ду 
TX (0,1) = т —5 (0, 0). 


et 
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Solution to SAQ 11 


PDEISAQII 


» 
x 


Clearly the string is continuous at x = 0. Therefore one boundary condition is 


y (0.0 = y* (0,0) (t 2 0) 


This, of course, implies that 


ген Ayt 
(00 = 0.1) (1 > 0) 
and 
927 A2 yt 
E (0.0 = 5-0.) {>й 


Let the tangents to the string at х = 0 make angles 7 апа y* with the x-axis on the 
negative and positive sides of x = 0, respectively. Consider a small element of string 


near the junction between x = —3Ax and x = $Ax. It has mass 
Мот sec ^ + р? sec *)Ax 


and acceleration 


y eyt | 
PL (0. o|- rm о) 


The y-component of the force acting on the element is 
T(siny* — sin р) 
Newton's Second Law of Motion now gives us 


zi 


T(sin у? — sin ут) = (p^ sec" + p — Jax Sa (0, 2) 


Proceeding to the limit as Ax approaches zero, we obtain 


siny* = sin ý”. 


ie. 

tan y* = tan №7. : 
But 

tan y — 
so that 


E 


This is the second boundary condition at the junction. 


( > 0). 
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Solution to SAQ 12 
The equation is 


ry 0 <x<imr>0). 


We have the fixed boundary condition 
(0,1) = 0 и 2 0), 


and the free boundary condition 


nn =0 (20). 
Cx 


The functions 

fixesSsinx — (xeR) 
and 

g:x 0 (xe А) 


are odd about x = 0 and even about x = 41, as required by the boundary conditions. 
Using the general formula for thé solution 


xte 


u(x. t) = M + ct) + f(x — cd) + x] g(X)dx. 


x-a 
we obtain, since c = 1, 
u(x, t) = $[sin(x + t) + sin(x — 1)] 


=sinxcost (tf 20,0<x < $n). 


Solution to SAQ 13 
The functions 

f:xesininx (0<x <1) 
and 

gixe—sinizxx (0<x <1) 


are odd about x = 0 and even about x = 1 and are therefore appropriate extensions 
of the initial conditions for the given boundary conditions. Using the general formula 
for the solution we obtain 


u = f(x + ct) + f(x — c0) + x[. d g()dx 


x-a 


2 
= sinizx cos 4nct + gn 3nx sin nct O<x<t,r20). 
лс 


Solution to SAQ 14 

There is a free boundary condition at x = 0 so that we require 
Дх) = /(—х) and g(x) = g(—-x) 

There is a fixed boundary condition at x = ! so that we require 
S(l+x) = —fü-— x) and g(l4 x) = —g(l — x). 

Putting — y = 1 + x in the second equation for f we obtain 
f(-») = -fQl + y). 

But 
fi) = fi-». 
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and therefore 
JO) = —f0l + у). 
Replacing y by (21 + y) in the last equation, we obtain 
fl y) = —f(2l + Ql + y) = —/(41 + у). 
The last two results give us 
J(y) = f(4l + y). 
Similarly, 
80) = g(4l + y). 


Thus, when we have a free boundary condition at one end and a fixed boundary 
condition at the other, the period of the extended initial condition functions is 4l. 


When both boundary conditions are free or both fixed the period is 2/. 


Solution to SAQ 15 


The only true statement is C. 


Solution to SAQ 16 


We can extend f and g as even functions about x = 1 and odd functions about x = 0 
by using the definitions 


f(x) = x? — 3x (-I€x&1) 
(2 = xP – 3(2-— х) (1€ x <2) 
=(2 + xP +3(22+x) (-2&x&-I1) 
f(x — 4n) (4n — 2 <x < 4n + 2, n an integer), 
and 
g(x) = [sin xx| (0<x < 2) 
—|sin zx| (-2<x <0) 


g(x — 4n) (An — 2 < x € 4n + 2, n an integer). 
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The Principle of Superposition 
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Introduction 

Second-Order Equations with Constant Coefficients 
Characteristic Curves 

Numerical Method of Characteristics 


Summary 


Solutions to Self-Assessment Questions 


РРЕ 2 


Set Books 


G. D. Smith, Numerical Solution of Partial Differential Equations (Oxford, 1971). 
H. F. Weinberger, A First Course in Partial Differential Equations (Blaisdell, 1965). 


It is essential to have these books: the course is based on them and will not make 
sense without them. They are referred to in the text as $ and W respectively. 


Unit 2 is based on. W: Chapter II, Sections 6 to 9. 


Conventions 


Before working through this text make sure you have read A Guide to the Course: 
Partial Differential Equations of Applied Mathematics. References to Open University 
courses in mathematics take the form: 


Unit M100 13, Integration II. for the Mathematics Foundation Course, 
Unit M201 23, The Wave Equation for the Linear Mathematics Course. 


PDE2.0 


2.0 INTRODUCTION 


In this text we are concerned with linear partial differential equations which have 
linear subsidiary conditions, and, in particular, we examine second-order equations 
in two variables. 


The text divides into three main parts. Section 2.1 concerns the linearity properties 
of equations and their subsidiary conditions, and introduces the principle of super- 
position for nonhomogenous problems. In the second section, second-order equations 
are classified into three types and characteristic curves are discussed. Finally, the 
numerical method of characteristics is used to obtain numerical solutions of hyperbolic 
equations. 


PDE 2.1.1 
2.1 LINEARITY AND SUPERPOSITION 


2.1.1 The Linear Partial Differential Equation 


In Unit 1, The Wave Equation we obtained the solution to the wave equation 


in the form 
u = p(t) + 401). 


where p and q are any twice differentiable functions of onc variable and € = x + ct, 
y = x — ct (see W: page 9). i 


In this section we are going to generalize this form and demonstrate that it represents 
solutions ofa wider class of linear partial differential equations. 


Suppose that a function (х, t) — u is given by 
u = р(&) + 40) + č, 


where p and q are arbitrary functions of one variable, and ё, 7 and ¢ are variables 
dependent on x and г, but not necessarily according to 


€é=x+e,y=x—ct,6=0 
as in Unit 1. 
The only constraint we impose on the prescribed functions 
(0 э 6, 
(x,t) —> 1, 
(x, J — C 
is that they be twice differentiable. 


We partially differentiate the expression for u, with respect to x to obtain ĉu/ĉx, and 
with respect to t to obtain @u/ér. By the chain rule for functions of several variables, 


Qu й дё [4 
Fr pm. + @ [rus D. x 
and 
2 on... at 
=p ой z qr + ap 


where p' and q' are the derived functions of p and q respectively. The second-order 
derivatives of u are given by 


д?и 1 дЕ} 022 „уд. 1026 
Е 6 + p) ates xj + (2 + o 


дх? дх 

^u „уб 0€ E » „дп ON д? 2; 

дхйг lay ae + Pss; 1 y, t TO Rae + axa 
апа 

ди, e „nE [On e 

95 = р"(&) Ё + Psa +4 (n) 2j + qi no з 1% am 
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We now havea set of five equations involving the four quantities p"(é), P'(£), q" (ij) and 
g'(n). Moreover, the equations are linear in these quantities, and may therefore be 
written in the form 


ec en OC ди ] - 4 
WE LE WE a per unt 
дё On 24 ди 
0 [i4 en Qt Ou - 

Gt д! й д pa : 
og? aE Ф\2 22у PE ди „ -lo 
ôx Ox? дх 6x? ox? ax? #0) |= 
0202 HE nay у ex д?и " 0 
dxdt бхй Ox O0 xð дхй Ox e 
4r fo (y omy C ош | | | | 
a ón — er o! а? а? | > | 


Since the column vector (p"(é), p'(£). q”(n)} 90). D* is not the zero vector, the matrix 
in the equation must be singular, and so its determinant is equal to 0. Now, č, n and € 
are known variables, and derivatives of u only occur in the right-hand column of the 
matrix. Asa result, when the determinant is expanded and set equal to zero. we oblain 
2, u 4 1 


O*u a? o ди ди 
at B(x, UA; + C(x, az + Dix, N= + Ext + F(x, 1) = 0, 


A(x, t) 
where 4, B, C, D, E and F are known functions of two variables. Thus 
u = pd) + qi + 


is the solution of a second-order linear partial differential equation and we can always 
obtain such an equation from an expression of this kind. 

It is important to realize at this point that we are not proposing a solution of the above 
form for all linear partial differential equations. We shall see in Section 2.2 that the 
equations of a particular, well-defined class have solutions of the above form. 


You should check that if £, у and č are given by 


(xccl 
n-x-—ct, 
¢=0, 


then we obtain the equation 


as expected. 


SAQ I 
If u = p(x? — t) + q(x? + 1). where p and q are arbitrary functions of one variable, 
prove that 


(Solution on p. 28.) 


* For typographical convenience an n-tuple may denote a n x | column matrix. 
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SAQ 2 
Show that, if p and q are arbitrary functions of a single complex variable. then 
u= pix + at + ify) + gtx — 20 — ify) 


is a solution of 


Ow ёи € 


ecco 7 
provided that 22 + ff? = 1. 


(Solution on p. 28.) 


SAQ 3 
Prove that 
u = [pt) + qu)! ^. 


where p and q are arbitrary twice differentiable functions of a single variable and 
where č, are variables dependent on x and t, satisfies a second-order (nonlinear) 
partial differential equation. 


(Solution on p. 28.) 


We shall now proceed to consider the linearity aspect in more detail. Before reading 
the next passage, however, you may find it helpful to review Unit M201 23, The Wave 
Equation, Section 23.1.2. 


READ W: Section 6, page 29 to page 32, line 21 (the end of the example). 


You may ignore the references in the two examples to Section 5 of W which is not 
included in this course. 


Notes 


(i) W: page 29, line —5 
Not only do most physical problems involve nonlinear operators the resulting 
equations are usually either very difficult or impossible to solve analytically. 
However, where they can be approximated by linear operators, solutions may 
well be obtained which give a good idea of the ways in which to tackle the non- 
linear problem. 


(i) W: page 31, line 6 
The solution set of every homogeneous linear partial differential equation is an 
infinite-dimensional vector space (Unit M201 23). Similarly the set of functions 
satisfying the homogeneous lincar subsidiary condition 


Liu] = 0 


is an infinite-dimensional vector space. in general. The general solution of the 
homogenous linear problem is the intersection of all these spaces, which may be 
finite-dimensional. 


(iil) W: page 31, Example 
The problem may be interpreted physically in terms of a string whose ends are 
fixed at x = 0 and x = 1, which is subject to a transverse force which varies 
along its length as sin zx. (Note that the ends are subject to no body force.) At 
time = О the string is released from a stationary position along the x-axis. 


The steady-state (or equilibrium) solution of the partial differential equation 
with the boundary conditions (but not the initial conditions) is given by r(x. 1) 
subject to 
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The equation now reduces to an ordinary differential equation whose solution 
(subject to the boundary conditions) is 


u(x,t) = — sin zx. 


If we put w = u — v, the problem for w may be interpreted as a string with no 
external force released from an initial displacement distribution 


io 
— —5 SIN zx. 
л? 
This problem is the vibrating string problem with c = l and / = 1, and its solu- 
tion is given by Equation (2.16) on W: page 13. When we superpose the steady 
displacement 
Ya 
u(x,t) = — sin zx 
E 
we obtain the standing-wave solution 
T 
uix. t) = -3 Sin zx(1 — cos nt). 
л 
This is readily seen to bea standing wave by observing that x = O and x = 1 аге 
points of zero displacement for all time. In other words, the wave has the form 


A(t) sin xx, where the profile does not change its basic shape with time — just 
its amplitude varies. 


5404 
If 
a2 a 
ENT Б agn = 2cosx — x sin x О<х<лп, г> 0, 
и(х, 0) = 0 О<х< л, 
10.0) =0 0<х<л, 
e 
(0.0) = 0 120, 
ulz, t) = 0 120, 


determine ($7, іл). 


(Solution on p. 29.) 


SAQ 5 


Which of the following are linear operators? 


E 


д a2 a 
Ou | Cu 1ди 
x АЗ — + cosu 
Ct ex? x Ex 


Lite 


(Solution on p. 30.) 
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SAQ 6 
If 
ё?и Ou 
EE (x,t) 52070 20 
ц(х,0) = 0 
ди 
ac 0-20 


find u(3, л). 


u(0, t) — sin? t 


ull, t) 2 0 


(Solution on p. 30.) 


x) cos 2t 


2.1.2 The Principle of Superposition 


PDE 2.1.1/2.1.2 


READ W: page 32, line 22 to page 35 (the end of Section 6), omitting the final sentence 


on page 32. 


Notes 


ü) W: page 33, line — 10 and footnote 
A series of functions converges uniformly if it is convergent in the uniform norm 
(Unit M201 19, Least-Squares Approximation). If a series of functions converges 
pointwise to g, a sufficient condition for term-by-term differentiation to be valid, 
i.e. for the differentiated series to converge to g’, is that the differentiated series 
converge uniformly. The footnote provides an alternative criterion for uniform 
convergence. 


Do not worry if this is not clear at this stage—it will be discussed again in 
Unit 6, Fourier Series. 


(ii) 


W: pages 33 and 34, Example 


Do not confuse i here with ,/ — 1! 


In this example v? satisfies the equations at the top of W: page 33, with the 
operators given by 


and 


The proble 


t 


L, 


L, 


\ 
лв 


а? ax?’ 
и] = х u(x, 0), 
al] = x (0). 
alu) = хэ u(0, t), 
и) = x utl, t). 
v(x) = sinimx  ieZ*. 
m specifies (W: page 34) 
u(x) = У, i7 * sin inx 


i7 * L po? (x). 
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Therefore, 


u= Yit 

і= 1 
is a solution, since the series obtained by applying the operators L, Li, Lz, Ly 
and L, to u converge uniformly—-a fact which we do not propose to prove here! 

(iii) W: page 34, line 21 

The uniqueness problem for the solution of the linear problem (6.1) on W: 
page 30 depends only upon the linear operators L, L,,..., L, and is independent 
of the particular data F, f,,... f. In other words, it depends upon the linear 
partial differential equation having specific types of linear boundary conditions, 
the actual functional forms of which do not affect the uniqueness of the solution 
if one exists. Uniqueness is determined without precise knowledge of the prob- 
lem sit is sufficient to know that the corresponding homogeneous problem has a 
unique solution. 


5407 
Show that the problem 
Ou T 20?и | Еі | 
ago = aoc 0 = F(x) 0<х<и1>0, 
u(x, 0) = f(x) O<x<i, 
д 
—wO-g) OS x <I, 
ôt 
u(0,1) = 0 120, 
д 
410) =0 120, 
ôx 


has, at most, one solution. 


(Solution on p. 31.) 


2.1.3 Uniqueness for the Vibrating String Problem 


We have seen in Unit / that the wave equation with constant coefficients can be 
solved, uniquely, for fixed and free boundary conditions. In terms of the initial con- 
dition functions for u and ди/д we found that 

xc 


u(x,t) = i f(x + ct) + f(x — ce) + al g(X)dx. 


xm el 


Clearly a change in faffects u(x, г) only atx — et or x + ct. Thus the effect ofan initial 
displacement is propagated at precisely the speed c. However u(x, t) is affected by a 
change of g anywhere in the interval [x — ct, x + ct]. Thus the effect ofa change in the 
initial velocity is propagated at all speeds up to c. 

Now consider the point (x,,f,) and the two straight lines 


X — ct = constant 
and 

X + ct = constant 
which pass through it. These lines represent propagation with the maximum speed c 
in the positive and negative x-directions, and are called characteristics of the wave 
equation. A pair of characteristics passes through each point (x,,¢,) in the domain 
of the equation. 
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© 
a 
Ж 
S 
Md 
KY 
SX 
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Ў 
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i 
“se 
5% 
д 
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ee 
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OQ? 
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LEK 
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0, 
C? 
KX 
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ks 
e? 
М, 


(Xa. 6) 


(xi, 6) 


“y 


If a disturbance is concentrated at (x, г,), it can influence only that sector of the half- 
plane t 2 t, formed by the two lines of slope + 1/c passing through the point (x,, t4). 
Hence this sector is called the domain of influence of (x, !). Similarly we can form the 
domain of influence of the point (x;, t2). Note that all domains of influence intersect 
for t > 0. 


The domain of dependence of the point (X, 7) is the set of all points with / > 0 whose 
domains of influence include the point (X, t). Thus the domain of dependence of a 
point is the domain cut off from the upper half-plane by the two downward drawn 
lines of slope + 1/с. Clearly, since the domain of the wave equation does not include 
t < 0, it is possible to have non-intersecting domains of dependence. If the domains 
of dependence of, say, (X, , /,) and (X5. 72) do not intersect, the displacements at these 
points will not be coherent —they will be caused by initial positions and velocities 
which are independent of each other. 


Gi) 


Gai) 


"Y 


Because the pair of characteristics through a point bound its domain of influence, it 
turns out that if there is a discontinuity in a derivative of the solution across a given 
characteristic at some point, then this discontinuity propagates along the characteristic. 
We shall not prove this result, or even formulate it precisely, but we shall provide an 
illustration. 
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Example 


Consider the problem 


a2 22 
Cu 26 и 
эз 025 = XER x the upper half-plane), 
à ad 0 (x. t)e К x R* (th half-plane), 


u(x, 0) = f(x) XER, 
Qu 
3 9 =0 xeR. 

As explained on W: page 10. the solution to this problem is given by 
u(x, t) = [/(х + ct) + f(x e) 


at all points in the domain for which f is twice differentiable. Suppose now that 


X x20 
fF:x— 
х? х< 0. 
Then 
х® + се? 0 € r « x/c 
u(x, t) 9 4 2erx t 2 |x/e| 


—(x? + с202) 0< г< —х/с. 


The functions u, би/ёх and ĉu/ĉt are continuous throughout the domain R x R*. 
However, the partial derived functions 


are not defined along the lines x — ct = 0 and x + ct = 0, their values undergoing 
jump discontinuities across these lines. 


These lines are precisely the characteristics through (0,0), the point at which the 
initial condition function f does not have a second derivative. We say that the dis- 
continuity (in f”) propagates along the characteristics. 


We turn now to the more general vibrating string problem in which the coefficient 
c?(x) is no longer fixed along the string. The next reading passage develops the dis- 
cussion of characteristics for this case in the context of proving that its solution is 
unique. 


READ W: Section 7, pages 36 to 40. 


Notes 


(i) W: page 36, lines 11 to —6 
This passage derives an integral form of the differential equation. Its physical 
interpretation is an application of the Law of Conservation of Energy to the 
string. This law is discussed in Section 8.1 of Unit MST 282 8, Work and Energy 11. 
Briefly, it states that the change in the TOTAL ENERGY of a system is equal to the 
work done on it by all the forces. d 


We can obtain formulas for the kinetic and potential energies by studying the 
element of the string at x with undisturbed length Ax, at time r. 


(ii) 


PDE 2.1.3 


The Kinetic energy of a particle is defined as 
J x mass x (velocity). 


For our element of string we obtain, approximately. 


Й ди j 
зро) a(x. | Ах. 
ct 


The potential energy of a stretched string (i.c. the energy stored in the string) is 
equal to 


lension x extension. 
The extension of our clement of string is, approximately, 

Ах)? + (Àu?]*? — Ax. 
Since we are considering the situation at a given time ( we find that the potential 
energy is, approximately, 


1/2 A 2 
ne ijs = точь r) Ax. 
ex 


on expanding the bracket by the binomial theorem and ignoring higher powers 
of @u/Ox, which is assumed to be small, as in the derivation of the differential 
equation in Unit 1. The TOTAL ENERGY of the whole string at time г is then 


In addition, work is done in moving the ends of the string. The rate of work 
being done is 


force x velocity. 


Ат cach end of the string the rate of working is therefore the y-component of the 
tension multiplied by the transverse velocity. 


The left-hand side of Equation (7.3) represents the total energy difference for the 
whole string at times 7 and 0. The Law of Conservation of Energy requires that 
this be equal to the work done on the system between these times by the external 
force. 


F gives the external force per unit mass, so F(x, r)p(x)Ax is the force on our element 
of string. Thus the rate of doing work on the clement is 


ди 
F(x, Deos Ax 


and, for the whole string, this yields 


| Lm | 
Fpdx|dt 
0 0 e 


W: page 37, line 14 
From 


Qv 
" = < 
a0 = 0 0< <i, 


we deduce that v(x, 11) = ro(x). Then the condition 
v(x,0)2 0 implies рох) = 0, 


since г к—э v(x, t) is continuous at 0. 
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(ii) W page 37. line —8 
Note that c depends upon x here. and we do not have characteristics which are 
straight lines as before. We therefore seek an alternative method of determining 
the domain of dependence of the point (х, f). 

(iv) W: page 37, line — 1 
Using the formula for integration by substitution (backwards method), the last 
two terms of the left-hand side become 


T æ д i дь дь 

-Í T ix Dx, t) dt +f To x. D(x, t)dt 
ta OX et n. Ox et 

(x.t)on C? Ix.r)on Ci 


* OB ‚уде " 
= -[ ToS Qalx) (x, 0,(х))0(х) dx 


П 


* Í тух, ADNE (x, О(х))О\(х14х 
о Üx ct 


1 д e n d 
= f тух.) (хи) ах. 
o 0х er dx 
tx.r)on C4 +C 
(№) W: page 38, lines 12 to 19 
Note that 


ILICE А Jayy -— 
КЕ ICT z] + |2 1а) Na] [^ 


C+C: 


for all C, and С›---іп particular, when C,, C; are given by 


NN 
ах) сх) 


In this case, the equation becomes 


[viz ^ eV ах + fi 

у — qe ах 3. 
IPRC CX HAT 
: e 


€ 


This is possible only if 


e e 
== +е;- =0 oC, 
д ax 
and 
e дь 
=- сз =0 onc, 
t ex 
SAQS 


W: page 40, Exercise 3 
(Solution on p. 31.) 
5409 

W: page 40. Exercise 4 


(Solution on p. 31.) 


PDE22022.1 
22 CLASSIFICATION OF SECOND-ORDER EQUATIONS 


2.2.0 Introduction 


We have seen that the solution of linear partial dillerential equations poses difficulties 
far greater than the solution of ordinary lincar differential equations. However, for a 
particular problem - (ће wave equation-— we were able to find a change of coordinates 
which allows us to obtain the solution in terms of initial and boundary conditions. 
In addition we were able to determine a uniqueness theorem for the "generalized" 
form of the wave equation in which the coefficient c is not constant. Our success in 
both cases exploits the existence of characteristics for the equations. 


We now procced to extend our approach to an arbitrary linéar second-order partial 
differential equation in two variables. 


Before continuing with this section, you may find it helpful to review Section 32.3 of 
Unit M201 32. The Heat Conduction Equation, which discusses how the plethora of 
second-order linear partial differential equations can be sorted into three categories, 
and the type of subsidiary conditions appropriate to each. 


2.2.1 Second-Order Equations with Constant Coefficients 
READ W: Section 8, pages 41 to 43. 


Notes 
(i) W: page 41, line —3 
We specifically exclude x = fj = Qand y = ô = 0,sothat,if A # 0, the equations 
on W: page 41, lines —7 and —6 imply that « # O and y 0. 
(ii) W: page 42, line 3 
This means that č and y should be linearly independent. 
(iii) W: page 42, Equation (8.5) 
In obtaining this equation. we required A # 0. If, in addition, C = 0. Equations 
(8.4) become 
ё = 2Ax, 


n = 2Ax — 2Bt, 


or alternatively, 


=x, 


which gives the same ratios B/x and б/у. The hyperbolic condition becomes 
B? > 0, which holds for all (real) B. 


The case A = 0, C # 0 can be deduced immediately from the symmetry of 
Equation (8.1), and we can write 


ё = 2Ct, 
n = 2Ct — 2Bx, 
or 
б =1, 
=x- S 
B 
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(iv) W: page 42, line —8 
We require 


Af? + ВоВ + Ca? = 0, 


or 
2 
АЁ, + p? +C=0. 
a a 
Since B? = 4AC, 
E В B 
A- + B-+—= 
a? a + 4A 9 
and 
в. в\? 
вара 
: Qu, 
The coefficient of is now 


2 


2 2 
ABS + В(об + fy) + 2 ay «|2402 + ala + 2) 4B } 
a 


A 
B В? 
= .{-вю + als - z) + zh 


=0. 


(v) W: page 42, line — 1 
The standard form of a linear partial differential operator (or equation) is more 
usually called the canonical form. 

(vi) W: page 43, line 2 
A = 0 = B = 0, since В? = 4AC in the parabolic case. 


W: page 43, Equations (8.8) and (8.9) 
Since we cannot make the coefficients of 2^u/0Z? or д?и/дп? vanish, we make the 
coefficient of 6?u/@é@n vanish instead : 


2ABS + Blað + By) + 2Cay = 0 


(vii 


or 
и у «у 
2А s+ = = 0. 
+ aff +3) ens 
We choose 
3 а 2A 
bm -5 


Thus (8.8) is the desired transformation, the denominator being chosen as 


J4AC — B? so that 5 


Qu au 
Ци] = A| туз + 55 |. 
eC en 


(viii) W: page 43, line 19 ' 
Laplace's equation has no (real) characteristics along which discontinuities 
would propagate. The conclusion drawn by Weinberger must be taken on trust. 


The general solution of Laplace's equation is not as straightforward as the others 
and will be studied in later units. 

SAQ 10 

W: page 44, Exercise 2 


(Solution on p. 32.) 
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SAQ 11 

W: page 43, Exercise 1 
(Solution on p. 32.) 
SAQ 12 

Reduce the equation 


Ou eu 
+4 


en? Oxet 
to standard form. and hence solve it. 


(Solution on p. 33.) 


2.2.2 Characteristic Curves 


Before asking you to read W: Section 9 on general second-order operators, we sum- 
marize what we already know about characteristics. Firstly, we know that they define 
domains of influence and dependence: secondly, discontinuities propagate along 
them. In order to look more closely into their nature, let us use the second-order 
operator 


Cu Ou д?и ди ди 
Liur A(x, N55 + Bix. t) —— + Сх, Daz + Flu. =} 
дг? Oxet Ox at’ Ox, 
We wish to consider the problem of determining the solution of the equation [и] = 0 
which takes prescribed values on a given curve Г in its domain. By a curve in R? we 
mean a subset of R? given by 


(x, 02x = (s), t = 005), se Г) 


where / is an interval on the real line and д, 0 are functions / — R. We usually 
restrict ourselves to a piecewise continuously differentiable curve: that is, one for which 
ф and 0 are piecewise continuously differentiable on J and there are only a finite 
number of points in / at which the derivatives of g and 0 vanish simultaneously. 

We shall specify Cauchy data on Г, i.e. x, £, u, &u/€x and ujët are given for se Z. In 
practice we do not specify both би/бх and @u/@r at any point, because и. би/ёх and 
ёи! cannot all be assigned arbitrarily along the curve. To show this, we suppose u 


is known along the curve, so that du/ds may be determined. Applying the chain rule 
we obtain 


du Qudx дий 

ds Oxds ^ Ot ds? 
and we see that ди/дх and ди/дг are not independent since dx/ds = q'(s) and 
dt/ds = 0'(s) are determined by the defining equations of the curve. For example, in 
the Vibrating String Problem (Equations (7.1) on W: page 36) we specified du/dt 
along the initial curve ((x, 1):0 < x € l, 1 = 0! but did not specify explicitly Qu/Ox. 
Now, on the curve Г, put p = ди/дх and а = Gu/dt, so that, by the chain rule, 


dp Oudx диа 


ds ёх" йз — Otx ds 


dq _ ди dx ude 
ds dxdtds ` OC ds 
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Combining these with the second-order partial differential equation L[u] = 0, we 


have the following system of equations for 22/012, @u/éxét, @2u/dx? evaluated on the 
curve Г: 

д?и dq 

дг? ds 

д?и 4 

M| — |= ар К 

дхдї ds 

д?и 

Ad =F 

L Ox 


where M is the matrix 


dt dx 0 
ds ds 
dt dx 
@ 22 2 
ds ds |" 
A B С 


and А, В, C. Е denote the respective values of these functions. Thus, if det M # Owe 
can easily solve for the second derivatives 0?u/ét?, 0?u/@xét, д?и/дх? on the curve 
since х, t, u and the first derivatives are specified. 


By successively differentiating these equations it may be shown that higher derivatives 
of u of all orders are uniquely determined at each point on Г for which det M # 0. 
The values of the function u at neighbouring points can now be obtained by using 
Taylor's theorem fer functions of two variables. Thus we are able to conclude that 
the solution of L[u] = 0 is unique in the vicinity of a curve on which Cauchy data are 
given, provided det M # 0. 


SAQ 13 
Find det M and show that for elliptic equations it is nonzero. 
(Solution on p. 33.) 


Let us now consider what happens when det M = 0, so that M is singular and we 
cannot solve uniquely for the second derivatives on the curve on which the data are 
given. Then 


dx]? dx dt diy 
А] —-B——-4C|—| = 0 
ds ds ds ds 
If di/ds # Othroughout the curve then for each value of t there is a unique point on the 
curve. Thus the curve determines a function 

tex, 
and by the chain rule we have 

dx dxdt 

ds dt ds 
Thus 


dx)? dx 

A|—-|-—-B--Cz-0 

(а) di 

A curve satisfying this equation is called a characteristic (curve) of the partial differen- 
tial operator L. (In the special case of the vibrating string operator we obtain the 
result on W: page 38.) If we are given values of u and its first-order partial derivatives 
on a characteristic, our problem does not possess a unique solution at points off 
it, however close they may be. 
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We may solve for dx/dt to obtain 


which is a pair of ordinary differential equations. In a region with B? > 4AC (hyper- 
bolic), there are two families of characteristic curves, whilst B? — 4AC (parabolic) 
results in one family, and B? « 4AC (elliptic) yields no characteristics at all. In the 
case where the coefficients are constants we have seen these results in Section 22.1. 
Note that, since A, B and C are now functions the operator L may fall into different 
classes for different regions. 


For hyperbolic operators, i.e. for regions in which an operator has two families of 
characteristics, we shall employ the theory developed before SAQ 13 to derive a 
numerical method of solving the Cauchy problem, in Section 2.3. For parabolic and 
elliptic operators, however, the method fails. This may be explained by virtue of the 
inapplicability of Cauchy data to these operators since they do not yield well-posed 
problems (Unit M201 32). 


READ W: Section 9, pages 44 to 46, ignoring the reference to Section 82. 


Notes 


(1) W: page 45, lines 2 to 5 
The two families of characteristics are the curves č = constant and y = constant. 
Differentiating along a curve € = constant given, say, by 


{(x, Dix = ets) t = 0(8), s€ 1). 


we get 
дей _ 
xds ds 
and, provided «1/05 z 0 along the curve. 
_ 06/08 
dt^ ds|ds — | x 


from which Equation (9.3) follows. Equation (9.4) is derived similarly. 


(ii) W: page 45, line — 12 
Note that this is everywhere hyperbolic. 
(ili) W: page 45, line —6 
This is the solution of 
dx _ 1+? 
di =I+x 
such that x = a whens = 0. 
(iv) W: page 46, line 20 
These coefficients have the same sign as A because the polynomial 42? + BA + C 
has no (real) roots by virtue of the ellipticity condition B? — 44C < 0. 


SAQ 14 
W: page 47, Exercise 3 
(Solution on p. 33.) 


SAQ IS 
Classify the equation 
Cu Du P ĉu 
tT dx "x R* 
x ex (ere Rx К", 


and reduce it to standard form. 
Hence obtair, a general form for the solution. 


(Solution on p. 34.) 
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PDE 2.3 
2.3 NUMERICAL METHOD OF CHARACTERISTICS 


In our discussion of characteristic curves of second-order equations we saw that, if 
u, би/дх and ди/д were specified on a curve, then we could evaluate the solution of 
the hyperbolic differential equation 


e a2 
и Cru 
"i + B(x, Na 


nen Oxét 


=0 (1) 


Ou Qu du 
+ C(x, 05-5 + F(x, t, u, =, — 
p at” Ox 
at points off the curve (provided that it is not a characteristic) by using Taylor's 
theorem. The object of this section is to utilize this fact and producea usable numerical 


scheme for solving hyperbolic equations with Cauchy data. 


д 
д 


The first thing to note is that we cannot, in general, evaluate all the terms ofan infinite 
Taylor series in a finite time! In fact to evaluate more than two is quite tedious; we 
choose, therefore, the simple approach of using the first-order Taylor approximation 
(Unit M201 14, Bilinear and Quadratic Fornis) 


ц 


(х+ Ах, ERAT) 


We must remember that Equation (2) is only an approaimation, and to be of any 
practical use the quantities Ax and At must be small. Applying Equation (2) to points 
along the given initial curve we can evaluate the solution at points off the curve but 
close to it. If we could also find p = ди/дх and q = du/ét at these points, we would 
then be in a position to repeat the whole process and generate the solution of the 
differential equation at points further and further away from the initial curve. A 
large part of this section deals with finding a way to evaluate p and qat points off thc 
initial curve. 


In Section 2.2.2 we derived the system of equations 
dt dx 0 Oy dq 
ds ds д? ds 
dt dx д?и ар 
0 o |= (3 
ds ds дхд! 4з ! 
e? 
A B CI s =F 
ox 


for the solution of Equation (1) on the curve !(х. г): x = ф()1 = Os). se Ij. 
The curve is a characteristic if the determinant of the matrix vanishes at each point, 
ie. if 
a2 dx Й 
Ах [2] — во, + cix 0. (4) 
di dt 


In this case Equation (3) does not have a solution unless the rank of the augmented 
matrix 


dt dx dq 
lo gy 1 
ds ds ds 
dt dx dp i 
0 d dod 
A B С -F 


equals the rank ofthe matrix M of Equation (3) (sce Unit M /00 26, Linear Algebra 111). 
Now det M = 0 means that 


rank (M) « 3, 
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so that any three columns are linearly dependent. In particular. this is true of the first, 
third and fourth columns, and so 


dt, dq 
d ° d 
det dx dp |=0. 
0 rst n 
ds ds 
A С =F 


Hence, we can say that a solution of Equation (3) exists along the characteristics given 
by Equation (4) only if 


dx dq , dp di P рх dt 


ds ds ` — ds ds ds ds — 
which, if dt/ds # 0 along the curve, is equivalent to 


dx dq dp „dx | 


^ dt ds ds * ds — » D 


Now for hyperbolic equations, Equation (4) gives two distinct real values for dx/dt. 
Thus Equation (3) gives rise to two equations from which we might hope to obtain 
dp/ds and dg/ds, which can then be approximated by finite differences to yield p and q. 


Let us write the solutions of Equation (4) as 


dx е dx 

— = f(x.t) and -— = g(x.t). (6) 
qon а 47s 

It will be convenient to speak of the /- or g-characteristic whenever we wish to 
discriminate between the two characteristics passing through a given point. 


For the remaining part of the argument we shall refer to the diagram, in which Г 
is a non-characteristic curve along which initial values for u, p,q are known. P and Q 
are points оп Г, with coordinates (xp, fp) and (xo. fg), which are close together, and the 
[characteristic through P intersects the g-characteristic through Q at the point R 
with coordinates (xp, tg). 


I 


If Equations (6) cannot be solved analytically— although we should always take 
advantage of reasonable analytical solutions when available—we can approximate 
them by using the first-order Taylor approximation once again, giving 


Xp — Xp = (ig — lp) Sp 
and (7) 


Xp — Хо = (In — logo 


ко 
[5] 
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where Je = (хь) and gg = g(xo.tg. Similar approximations to dp/ds and 
dq/ds yield the two approximations to Equation (5): 


(dr — qp) (Pr — Pp) (xp — xp) 
Ap fr +C E 22 
i "sa = sp) "(ss — sp) Н я — 5р) 0 
and 
ldr — do) (Pr — Po) (Xr — хо) 
обо; — Sq) MR — 50) tra (Sr — 50) 


where again the subscripts P, Q and R indicate the values of the respective functions 
and variables at P, Q and R. 


Multiplication by (sk — sp) and {Sr — Sg) respectively gives 

Ap feld — dp) + Съра — рь) + Fp(xg — Xp) = 0 (8a) 
along the f-characteristic, and 

Ag8aldr — йо) + Сора — ро) + Fo(xg — xg) = 0 : (8b) 
along the g-characteristic. 


Since p and q are assumed known along Г and hence at P and О, and since we can 
calculate хь from Equations (7), Equations (8) are two simultaneous algebraic 
equations for рк and qg. Hence we have a method for calculating p and q at points 
off the initial curve and also, using Equation (2), to calculate и there. 


This process can be applied to the points Q and W to give the values of u, p and q 
at the point S. Hence, knowing u, p and q at both R and S, we can apply the procedure 
once again to evaluate u, p and q at T. Therefore, by taking more points along the 
initial curve we can progressively estimate the solution of the differential equation at 
points further away from the initial curve. The points off the initial curve are given 
by the intersections of the characteristic curves which pass through the chosen 
points on the initial curve. 


Note that in this scheme it does not matter whether we calculate p and q before or 
after we calculate и at a particular point. However, if we calculate p and q first then 
we can use an approximation for u which is better than Equation (2). To see this, 
we first write down the first-order Taylor approximations to Ou/Ox and du/ér at 
(x + Ax. t + AD): 


a e 8? 
ди ĉu Ou eu 
mo + Ах) + At (9a) 
OX|saxisAn Хм CX" |o. n CXCL s n 
and 
A 2 
Qu ди eu 
Е =f Axa ч А5 (9b) 
C |х+ахл+м) C! jon EXOT a, n CU juin 


Multiplying Equation (9a) by Ax and Equation (9b) by Аг, and adding the results, 
we obtain 


] TA E 
(хм) e 


ёх 


үх+Ахи+м) СА QURAxU + м) 
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The left-hand side of this equation is just the expression for the second-order terms in 
the second-order Taylor approximation : 


А 52 
ди ди ди 
и =u Ах м—| + МАХ)? I 
(х+ Ах + an (n OX fixan Cl fn X jt n 
2* aoo 
+ AxA HASS 
Oxt |a. n OU sin 
ETS E ди | 
=u 3Ax| — - 
wn OX|esascsan ÊX fin 
Qu 


+ 4Ar É 


et N 

This apparently first-order approximation is more accurate than Equation (2) since 
it incorporates the second-order terms of the Taylor expansion. Applying this 
approximation to the points P, Q and R gives the approximations 


(asreag OF 


Ug = Up + (рк + po) xg — Xp) + (4а + др) (а — tp) 
and 
ик = и + Xp + Ро) (а — Хо) + ida + dollir = tQ). 


We can use either of the last two equations to calculate ug, the remaining equation 
being used as a check on the method. 


We summarize the method as follows. Given a second-order hyperbolic linear partial 
differential equation in и with (initial) values of u, ди/дх, Ou/Ot given on a non- 
characteristic curve, the solution и can be found at points off the initial curve by 
repeating the following calculations: 


Stage (i) 
calculate the coordinates of the point of intersection of two characteristic curves 
given by (xg, tr) from 


Xp — Xp = (Ig — tp)fp. 
Xp — Xo = (tn — tago 
where f (x, t) and g(x, t) are solutions of the equation 
dx|? dx 
лә) — B(x, Dac + C(x, = 0; 
Stage (ii) 
calculate pp and qR from 
Аъ — de) + Ce(Pa — Pp) + (ха — Xp) = 0 
and 
AQSaldr — 49) + Со(рв — Pa) + FolXr — Хо) = 0; 
Stage (iii) 
calculate ug from both 
Mg = Up + (Pr + рв) (Хв — Xp) + (qn + qp) (fa — tp) 
and 
Hg = Ug + Ура + ро) (хв — хо) + lar + qo) (fa — to). 


the second calculation being used as a check on the method. 


Example (S: page 117, Exercise 1) 
The equation 


au au a +1<0 ey 
ax й “are 5 SESE ovd 
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satisfies the initial conditions 


eu 
и = z = хајопвг = 0,0 < х <1. 


PDE 2.3 


Use the numerical method of characteristics to calculate the solution at the points 
R, S and T shown on the diagram. The coordinates of the initial points P, Q and W 


are (0.4,0), (0.5,0) and (0.6,0) respectively. 


^ 


Solution 
The differential equation 


e m д?и OU 
Ox? дхй Сә? 


has characteristics specified by the roots of the equation 


which are 


dx 
—-i-l 


dt 
The characteristics are given by 
2x — t — constant, 


t + x = constant. 


Stage (i) 
We calculate xg and rg from - 
2х8 — fg = 2Xp — fp. 
Ig + Xp = lo + Хо. 
Stage (ii) 
Py and qg are obtained from the equations 


— (dr — qp) + (ра — pp) (xg — Xp) = 0, 


2(4к — do) + (PR — po) + (Xe — Хо) = 0. 
The values of p, q at P and Q are obtained from 


ĉu Ana 
q(x, 0) = a0 =® O<x<il, 


m 
a 
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and 
p(x.0) = C x. 0) = duc. 0) along t = 0 
Ox dx 
=] О<х& 1 
Stage (iii) 


The formula(s) for ug can now be used since wc have already determined x, t, p. q 
at the points P, О and R. The values of u at P and Q are given by the condition 


u(x,0) 2 x O<x<i. 


We now repeat the whole process for the points S and T, obtaining the data 
set out in the following table. 


exact value 
E d t р q и of A 
— = — 
P 04 .00 10 04 04 0.4 
Q 0.5 0.0 1.0 0.5 0.5 0.5 
№ 0.6 0.0 1.0 0.6 0.6 0.6 
R 13 2 32 15, 418 418 
30 30 30 30 900 900 
16 2 32 18 514 514 
3 30 30 30 30 900 900 
T 14 4 34 18 484 484 
30 30 30 30 900 900 


The exact values of u have been determined from the analytical solution 
u(x t) 9 x + xt 02/2 

from which we obtain 
р(х.) = 1 + tand д(х.1) = х + 1. 


The precision of our numerical solution is due to the fact that the quadratic 
approximation is exact. 


SAQ 16 

Use the numerical method of characteristics to obtain u (0.15, 0.15), where 
ðu д?и 
me = 17 О<х< 1, 1> 0, 


subject to the boundary conditions 
(0,0) = u(1,1) = 0 120, 
and the initial conditions ` 
u(x, 0) = 3x(1 — x) O<x<l 


eh = 0xx«l 
a9 = &х<1. 


(Solution on р. 35.) 
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2.4 SUMMARY 


In this unit, we have looked at linear partial differential equations of the second order, 
subject to linear subsidiary conditions. We showed how their linear properties could 
be used to give us the principle of superposition, and also to examine the uniqueness 
of solutions to these problems. 

Second-order operators of the general form 


Ou д?и ĝu 


Qu ди 
AL —- — х = 
PE +В + Ср + Fen | 


дг ax 

are classified as hyperbolic (B? > 4AC), parabolic (В? = 4AC)and elliptic (B? < 4AC). 
We discovered that for hyperbolic equations, if x, г, u, and ди/дх or Qu/8t are specified 
on a given curve Г, our problem possesses a unique solution throughout the region 
bounded by Г and a pair of characteristics. 


We found that, along the characteristics, x and t are related by 


1х\? > 
al og с шо, 
dt dt 


and saw how the equations of the characteristics can be used to find a coordinate 
transformation reducing the differential operator to standard form. 


Finally, we saw the use of characteristics in obtaining numerical solutions of hyper- 
bolic equations. 
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2.5 SOLUTIONS TO SELF-ASSESSMENT QUESTIONS 


Solution to SAQ 1 
We are given 


и = p(x? = 1) + а(х? + 0). 


So, 
Qu es ЖР 
>= = 2х[р(х® — г) + q'ix? +9), 
дх 
д?и ЕР "E Me e 
52 = 4x?[p(x? — 1) + д" + 0] + 2D G7 — 0,4 (х + 0]. 
апа 
“ШЖ "(x? — 0 + 4"(х? +9) 
EE =p 1 ч . 
Непсе, 
Ou lôu „ёи _ 0 
ax? хх OP 


Alternatively, if you are feeling energetic, you can multipiy out the determinant of the 
matrix on p. 7 and set it equal to zero. 


Solution to SAQ 2 


If 
u = p(x + at + iBy) + q(x — at — Ву), 
then 
eu ns . >. | 
gc P erat iBy) + q"(x — at — ifly), 
Qu c | . | 
à = cB + at + ipy) + ax ог — 180), 
u a, | , | 
wc + at + iby) + g'(x — at — ify). 
Непсе, 
ди Ou . 
— +5 = (1 [р(х ву ЕРТЕ 
ag + ae (1 — В)" (х + at + iy) + q"(x — at — iflyy] 
_ Ou 
| 007 


provided 1 — f? = à? ога? + f? = 1. 
Solution to SAQ 3 
We have 


u? = p(é) + 4(п). 


We may differentiate once, using the chain rule, to obtain 


and 
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Further differentiation of these equations yields 


du\? д?и ac? are дп\? i 
=p" "E roln, graye 

(=) + 0—3 rof eri anfa +d. 

оди ди Ou n 06 0€ 02 


axar  axar = PRG Р?) 


д? 


„у дд 
+40) 211 xr 


дхдї дх д 


+ q'(n) 


and 


Ou? | Pu, fa "m p? д? 
X + Quay = р «(8 + P; «ez + dins. 


for all values of x and t such that the second partial derivatives exist. These five 
equations may be expressed in matrix form as 


Г де дп ди "Iu 
0 x 0 x ~ us | [p (0) [0 
дё ôn ди n 
0 á 0 9r -uz р) 0 

Qn? 022 Oo]? Pn Qu? Pu " 
[2 a dej xd 7 2] -uze | eee 
д5 д5 are д] ôn Ou оди ди 2 д?и ( 0 
dxdt дхй дхй дхй axo | | TO 

ce? 2 Op)? ah ди\? ди 

> ii =2 = 
AU дг? ді дг? BU 202 L | | d | 


As in. the text (Section 2.1.1), the existence of the nonzero column vector 
(p" (E), р'(@), q" G1) q' (m, 1) satisfying this equation means that the determinant of the 
matrix must be equal to zero. 


On expansion, this results in an equation of the form 
ди ди ди\? д?и 
А(х, ШЕ + B(x, Duz + C(x, ШЕ + | 


ди ди Oru Qu)? | ^u 
X, t) z— = — ot) i — | = 0, 
+ D(x, qx à + “| + Е(х (6) + бта 
where A. B, C, D and Е are known functions of two variables. This equation is поп- 
linear unless C, D and E are the zero functions, in which case it is first-order. 


Solution to SAQ 4 


A particular solution of the differential equation is the steady-state solution, which 
satisfies : 


v(x,t) = v(x,0) 0<х «л, 1 20 


and 
eu r 
— (x1) = 2cosx — xsin x 0cx«mt»0 
дх? 
Integration yields 
ov ; P 
——(x,t) = 2sin x + xcosx — sin x 
= sin x + x cos Х 0О<х<л,г>0 
апа 
ціх.) = —cosx + xsin x + cosx 
= xsinx erek г> 0. 
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We have chosen constants of integration such that v satisfies the boundary conditions 
on uH. 


Next we let w = u — v: then w must satisfy 


ae hw 
Dro O0cx«mt-0 
oO ёх? 

w(x.0) = xsin x 0<х<л, 
ор = 0 О&<х& л, 
(0, 0) = 0 120, 

w(x, t) = 0 120. 


By Equation (2.7) on W: page 10, we have 
w(x, t) = Mix + Dsin(x + t) + (x — t) sin(x — 0)) 


in the triangular region 


Єх, t<n-x, 1270. 


Now 
u(x, f) = w(x, t) — x sin x: 
hence 
тл л 
и|=,=| = –-= 
272 2 


In fact, this result is independent of the boundary conditions at x = 0 and x = л. 


Solution to SAQ 5 


C and D. 
А А г ди i Qu 
A is nonlinear owing to the presence of the term u—, B owing to à 


2 

* and E owin 
Ox 8 

to COS и. 


Solution to SAQ 6 


The point ($, $7) is outside the triangle in which the solution is determined by the 
initial conditions. However, we cannot apply the formula solution — Equation (2.16) 
on W: page 13— directly because the differential equation is nonhomogeneous and 
the boundary condition at x = 0 is neither fixed nor free. We proceed as in the 
example on W: page 32 and define the function v, which satisfies the last two con- 
ditions, by 


vix, t) 9 (1 — x)sin?t 0<х< 1,120. 
Putting w = u — v, we obtain 


Ow Pw u Pu 
= = жос ER 2 
17} WP a 2(1 — x)cos2t 


=0 О<х<1, t»0, 


w(x,0) = 0 0<х<1, 


м(0, г) = 0 120, 
w(I,r) =0 120, 
whose solution is clearly 


w(x,t) = 0 0<х 


^ 
v 
© 
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u(x=(l—x)sin?s Ox <1, 130. 


в. 
з 
a 
a 
Il 
tà 


Solution to SAQ 7 
Let u, and u, be two solutions and put v = и, = uy. Then 


22 2. 
др _ Eu D 


=0 0cx«Lt-0. 


v(x, 0) = $6.0) =0 о<х<1 


000, 1) = (Һи) = 0 120. 


The formula solution — Equation (2.16) on W: page 13— gives v(x,1) = 0 so that 
и, = uz. Thus there is at most one solution to the original problem. 


Solution to SAQ 8 


Equations (7.7) on W: page 38 give the characteristics in the form 


= | dé 
FEE 


and 


А * dé 
Carat f IFE 


where (x, i) = (4, 3). 
Thus, on C, we have 
123 — tan^! $ + tan^!x 


and on C; we have 


t-3-t12n^'i-tan^!x. 


Therefore the domain of dependence of (1, 3) is the region determined by 
Oxx«i Ox«r«3-tan'$ctan'!x 


and 


«x«l, O«t«3-tan'!i-tan^!x. . 


Solution to SAQ 9 
Applying the same method as that used in the:solution to SAQ 8 we get 


1 1 


Сүй 2tan^! x — ап! 
and 

Coit = an^! $ — tan! x. 
The domain of influence of (5. 0) is determined by the two upwards characteristics 
and is given by 


1 1 


O<x<4, rzin 'ji-tan'!'x 
and 

11 
A 


Lex<l, ¿> tan! x – tan” 
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Solution to SAQ 10 


We have 


B 
х= х= стапа г -r. 
2A 


Therefore 
ex _ 1 
ex 2А” 
er et 
— = 0,—— = 1. 
Ox et 


Differentiating u according to the chain rule for several variables, we obtain 


eu B ди ёи 
п 2A ax’ б” 
Ou ди » 
ôx Ax 


au В? ди В ди дш 


Pu В ди x u 
хіі — 2A0x? да” 


In terms of the new coordinates we obtain 


? y eu ди В? u Oru Qu 
Ци) = — == — Bz— + А5 - 75 Dt CI 
Ud = 213: — Pawar + Aap? ~ 2a ax? + Pawar + Сах? 
eu B? Ou 
= A- |—— – С|-—з. 
д? |2 Jes 
Since B? — 4AC > 0). this has the form of the wave operator, 
д?и „ё?и 
Ци) = [Ss -¢ Avi 
with 
› _ B!-4AC 
ez y~ 
4А? 
Solution to SAQ 11 
(a) B? — 44C = — 3. The operator is elliptic and has no characteristics. 


(b) В? — 4АС = 0. The operator is parabolic and the characteristic through 


(0, 1) is 
х-20+2 = 0. 


(c) B? — 4АС = 12. The operator is hyperbolic and the characteristics through 


(0, 1) are 
х+ Qc 34-12-20 
and 


x+- At- 1) = 0. 
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Solution to SAQ 12 


We have seen in the solution to SA 
on W: page 42 the transformation 


Q 11 that the operator is parabolic. By the result 
€ = 2x — qt 
n-t 
transforms the equation into the standard form 
д?и 
ap^? 
whose solution is given by 
u = p(5) + nq(). 
Thus 
u = px — 4t) + tq(2x — 41), 
where p and q are arbitrary functions of one variable. 
Solution to SAQ 13 
We obtain 
ixi? x » 
det M = А2) вай сја)" 
ds ds ds, 
For an elliptic equation 44C — B? > 0; hence А # 0. Therefore 
dx|? dxdt — B?[d|? — B?[dr? dt\? 
dtm = | CEN sas] - aa Ы c 
ах Ва]? В?\ jdt\? 
ka * (6-2) 
z 0, 
since 44C — В? > 0 and dx/ds, dt/ds are not simultaneously zero at any point. 
Solution to SAQ 14 


Using the method of W: page 45, the characteristics are determined by 


бех + t?) + d + 7)? — 42е) 
er + 1? + (ех — (0) 


The condition 
dx/dt = ех 


gives 


and 
t=¢ e 
Let x = a when t = 0: then 


a 


emet, 
and hence 
а = = In(e^* + t). 


We therefore take ' 


č = — In(e^ +0). 
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The condition dx/dt = 12 gives 
х= di c. 

Let x = a when r = 0: then 
fg = а, 


and hence 


We therefore take 


=x- 19. 


Solution to SAQ 15 " 
(2xt)? — 4x71? = 0, 


and so the equation is parabolic throughout its domain. Along a characteristic 


dx t 

d ox 
and 

x? = —1? + constant. 
So let 


ё= х? +, 
and take 
nt. 


We obtain 


o dë Om 
Qu ди 
ao Do 


Ou _ „ди 4, 2t Qu | Qu 


ar ag * 4 ар + Pp дү? 
Oru Uu д?и 
ахат 7 41902 + a 


The equation now becomes 


Qu д?и Ou ди д?и ди 
x12 4? 4 — 2) 2x—— 
x | at^ 222+ tmn + ET | ч[а + са 


which simplifies to 


u х? ди 


X =——. 
Oy t ôq 
Since x # 0 we may write this as 


A 1 ди _ 
д? n an 
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or 
[1 eu 0 
т\т On] — 
10и, 
" дп x: AG) 
or 
ди Т 
on nfi): 


лае OPA) + (8). 
Thus the general solution is 


u = Uf (х2 +P) + foe + 12). 


Solution to SAQ 16 
The characteristics for the operator are 


t + x = constant and г — x = constant. 


0 0.1 0.2 03 04 05 1 х 


Let P, О, W, V. К, S, X, T, Y, Z be as shown in the diagram. The lines shown are 
the characteristics. We calculate ug as follows. 


Stage (i) 
We calculate xg and tg from 
fg — Xp = fp — Xp, 


Ig + Xp = fg + Хо. 


Stage (ii) 
py and gp are obtained from 
dr — dp — (Pr — Pr) = 9, 
— (dr — do) — (Pr — Po) = 0. 
The values of p. q at P and Q are obtained from 


Qu 
qix.0) = 3 (x0) = 0 0<х<1 


апа 
А 
ріх.0) = EET) - B 0)  alongr =0 


= -х +3. 


ә 
л 


Stage (iii) 


PDE 2 SAQ 16 


ug may now be determined using the values obtained for хв, tg. Pr: 4r- up and 


Ng are given by the condition 


«d х= <209 


u(x, 0) = 3x(1 — x) йсх&1. 
We repeat the whole process for the other points, obtaining the following table. 
| analytical 
* ! P 4 К solution 
0 0 0.5 0 0 0 
0.1 0 0.4 0 0.045 0.045 Initial 
0.2 0 0.3 0 0.080* 0.080 values 
0.3 0 0.2 0 0.105 0.105 
0.05 0.05 0.45 —0.05 0.0225 0.022 
0.15 0.05 0.35 —0.05 0.0625 0.062 
0.25 0.05 0.25 —0.05 0.0925 0.092 
0.1 0.1 0.4 —0.1 0.040 0.040 
0.2 0.1 0.3 —0.1 0.075 0.075 
0.15 0.15 0.35 —0.15 0.0525 0.052 


7 


The analytical solution is 


u(x, t) = 


1s 


1 


1 т 
(1 — cos пл) cos пл! sin nx. 
n 


The first three terms of this series have been used for the comparison column in 


the table. 


In this question we have not had to use the given boundary conditions since we have 
calculated the solution at points where the solution depends solely on the initial data. 
However, if, in this question, we had wanted the solution at the point (0.05, 0.15) we 
would have needed to use the given boundary condition along x = 0. The chosen 
characteristics intersect each other at points on the boundaries x = 0 and x = 1 and 
therefore it is possible to "pick-up" these known values and use them in the numerical 
scheme. (If the chosen characteristics do not intersect at points on the boundaries 
we could still use the method of characteristics in the regions affected by the boundary 


values by choosing a new set of characteristics there.) 
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Set Books 


G. D. Smith, Numerical Solution of Partial Differential Equations (Oxford, 1971). 
H. F. Weinberger, A First Course in Partial Differential Equations (Blaisdell, 1965). 


It is essential to have these books; the course is based on them and will not make sense 
without them. They are referred to in the text as S and W respectively. 


Unit 3 is based on W: Chapter III, Sections 10 to 13. 


Conventions 


Before working through this text make sure you have read A Guide to the Course: 
Partial Differential Equations of Applied Mathematics. References to Open University 
courses in mathematics take the form: 


Unit M 100 13, Integration 11 for the Mathematics Foundation Course, 
Unit M201 23, The Wave Equation for the Linear Mathematics Course. 


Bibliography 


M. R. Spiegel. Vector Analysis (Schaum. 1966). The first two chapters contain some 
basic material about vectors. Chapter 4 discusses scalar and vector fields and intro- 
duces the gradient and divergence operators. Surface and volume integrals are 
covered in Chapter 5 and the Divergence Theorem is treated in Chapter 6. 


C. A. Coulson, Wares (Oliver and Boyd, 1955). You may find this book generally 
useful as it gives many examples of physical situations which give rise to the wave 
equation. The derivation of the membrane equation appears in Chapter III. 


PDE 3.0 
3.0 INTRODUCTION 


This unit is concerned with the general properties of two important partial differential 
equations which occur frequently in applications to physics and engineering. 


The first is Poisson's equation in two or three dimensions, 


a2 д2 
Ou @?и 
aoe Pag = F, 
ex? ду? 
or 
Qu Qu ey 
AT m3 3 F. 
Ox ey Oz 


Poisson's equation is the simplest example of an elliptic equation, and an important 
application of it is discussed in Section 3.2.2 for the flow of a viscous fluid in a pipe. 
The special case in which F is the zero function is Laplace's equation, which has-many 
applications including: 

(a) gravitational potential; 

b) steady heat flow; 


c) flows of non-viscous fluids; 


( 

( 

(d) electrostatics; 
(e) surface waves on liquids; and 
(f) elasticity. 


Uniqueness and maximum values of the solutions are discussed together with an 
explanation of what we mean by a properly posed problem. Solutions to Poisson's 
equation can be obtained analytically for simple boundaries and boundary con- 
ditions only. Section 3.2 introduces extremum principles, which can be used for 
obtaining numerical upper and lower bounds for certain integrals involving the solu- 
tion. Bounds for the flow of a viscous liquid in a pipe are presented in detail. 


The second partial differential equation which we discuss is the diffusion equation, 
which for three spatial dimensions has the form 
ĉu k ди Ou ul) _ 0 

ec ey gas] C 
It is an example of a parabolic equation, and one of its applications, which we model 
in this unit, is to the time-dependent conduction of heat. We discuss both the unique- 
ness of the solution and a maximum principle for this equation. 


We include some results in vector calculus where these are necessary for our develop- 
ment of the material; and some elementary facts about vectors are collected in the 
Appendix. 
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31 POISSON’S EQUATION 


3.1.1 The Membrane Equation 


In order to discuss problems in two and three spatial dimensions, it is convenient to 
employ vector notation; that is, we describe physical quantities (where appropriate) 
by geometric vectors, or by corresponding ordered pairs or triples. To remind you of 
the sort of vector algebra which we use we have set SAQ 1; refer to the Appendix if 
you are unsuccessful in interpreting it. 


SAQ 1 

Ifa =i+j+k,b=2i-j—kande = —i + 3j + k, find 
G) lal 

(ii) |a + 2b| 

(iii) a-b 

(iv) ax b 

(у) а:Ъхе 

(vi) a x (b x c). 

(Solution on p. 34.) 


We now proceed to the definitions of vector calculus which are required to enable the 
use of vector quantities in differential equations. In applications our domain of interest 
may be the plane R? or three-dimensional space R?. 


Let D denote a subset of R?, for example; then a function 
ф:р — R or o:D x R —>R 
(x, y, 2) эи (x, y, 2,0) и 
is called a scalar field. Thus a scalar field is just a generalization of the type of function 
encountered in ordinary differential equations, i.e. with domain and codomain a 
subset of the reals. We define addition of scalar fields pointwise over the domain as 


follows. If ф, y are two scalar fields with domain D c R? then ф + y is the scalar 
field defined by 


(6 + W(x, y, 2) = ф(х, y, 2) + WO, у, 2). 
Similarly we may define scalar multiplication by 
(тф)(х, y, 2) = т(ф(х, у, 2) — meR. 


With respect to these definitions the set of all continuous scalar fields with a given 
domain form a vector space. 


A vector field is a function v with domain D (in А?) or D x R$ and codomain G? (the 
vector space of geometric vectors in the plane). Thus 


v:G y) — vo yji + v(x, yj 


where v, and v, are scalar fields with the same domain as v. We can also define three- 
dimensional vector fields R? —> G?. 


By defining addition and scalar multiplication of vector fields as for scalar fields we 
again obtain a vector space of vector fields on a given domain. 


The temperature in a block of metal is an example of a scalar field, and the local 
velocity in a liquid determines a vector field. 


The gradient of a scalar field ф is the vector field 


дф. дф. ð 
gió cie Oj + Pk 
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The divergence of a vector field v = (vy, Vy, v.) is the scalar field 


Note that 


22 2? 32 
divgrad 6 = vig =h , 99 7 

ex^  80y* дт 
If S, V denote the vector spaces of scalar and vector fields respectively on some given 
domain, then grad:S — V and div:V — S are linear operators (vector space 
morphisms). We have written grad for the gradient operator to indicate that grad ф 
is a vector field; Weinberger, however, writes grad. We can define two-dimensional 
versions of grad, div and V? by omitting the z-dependence. You should note that é/ét 
does not figure in these operators even for time-dependent fields. 


Our definitions of grad and div presuppose a rectangular Cartesian coordinate system. 
Definitions which avoid this presupposition are possible (see Unit MST 282 7, Work 


and Energy I. Section 7.1). However the formal definition given above is sufficient 
for our purposes. 


The directional derivative of the scalar field at (x, y, 2) in the direction of the unit 
vector e is defined to be the real number 


grad $ ` e, 


which is often denoted by d¢/ée. It measures the rate of change of $ in the direction 
of e. Since e · grad ф = |grad ф| cos 0, where 0 is the angle between e and grad ф, we 
see that the greatest rate of change of $ occurs when e is parallel to grad ф. 


grad ф 


ф(х, у.с) = constant 


The solution set of the equation 

P(x, у, 2) =k б 
represents a surface in А? for each value of the constant К. If the unit vector e is 
tangential to the surface we expect е grad ф = 0, that is, the rate of change дф/де is 
zero since @ is constant on the surface. Since this is true for any tangent it follows that 


grad is normal (i.e. perpendicular) to the surface at the point. As we might expect, the 
greatest rate of change of ф occurs in the direction normal to the surface. 


For two-dimensional problems the solution set of 
ф(х, у) = К 
represents а curve іп R?, and the previous discussion is equally valid. 


In general we call any vector t such t · grad ф = 0 at the point P a tangent vector to 
the constant value surface (or curve) of ф passing through P. The set of all such t forms 
the tangent plane (or line) to the surface (or curve) at P. The unit normal vector at P 
to the surface (or curve) is clearly given by (grad ф)/егай | evaluated at P. 
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SAQ 2 


Find the gradients of the following scalar fields. 


(а) (х, у, 2) 9 (х2 yz 


(b) 
(c) 


(x,y, 2) 9 хус 


(x, y, 2) x? + yz + 23 


(Solution on p. 34.) 


SAQ 3 
If F:(x, y, z) (2xy + 23)i + x?j + 3xz?k, find div F. 


(Solution on p. 34.) 


SAQ 4 


i-j-k. 


(Solution on p. 34.) 


READ W: Section 10, page 48 to page 50, two lines below Equation (10.2). 


Notes 


(i) 


(ii) 


(iii) 


W: page 48, lines 1 to — 10 

The membrane may be thought of as a thin sheet of rubber stretched over a 
frame as in the case of a drum, or as a soap film on a metal frame (see cover). 
Ideally the material is perfectly flexible in that it does not resist bending but does 
resist stretcning. 


The membrane equation quoted on /ine — 10 is the natural analogue in two 
dimensions of the equation for the vibrating string derived in Unit 1, The Wave 
Equation. In fact, if c, F and и are independent of y, then this equation (line — 10) 
reduces to the one-dimensional vibrating string equation. The detailed derivation 
of the membrane equation is complicated and requires consideration of the 
curvature of a surface. An elementary treatment can be found on page 43 of 
Coulson, Waves. 


W: page 49, line 7 
The second condition is equivalent to 


ye) = (t9) 7 
x(t) x'(t) 


if x'(t,) # Oand х'(т,) з 0. This ensures that the curve C is smooth at its common 
point (x(to), y(xo)) = (x1). »(z))). 

W: page 49, footnote . 

For example the interior of the circle x? + y? = 1 isa simply connected domain. 
On the other hand the domain bounded by the circles x? + y? = 2 and 
x? + y! = 1 is not simply connected. A domain іп R is just an open interval. 


A domain is not necessarily bounded. For example, the subset {(x, y):x > 0, y > 0} 
of R? is a domain. 
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simply connected 


not simply connected 


(iv) W page 49, line — 12 to page 50, line 17 
A knowledge of electromagnetic theory is not required for this course. Ifyou find 


the next passage difficult in that respect, read through it and extract just the 
mathematics. 


(v) W: page 50, line 15 
In this two-dimensional case, the tangent vector 


lie in the same line. Thus t x grad u = 0 produces the required equation. 


General Comment 


A boundary value problem in which w satisfies Laplace's equation in D, and is 
specified on С. is known as a Dirichlet problem: if, alternatively, ĉu/ĉn is specified 
on C, then it is known as a Neumann problem. Sometimes a mixed problem may be 
given, in which u is specified on part of C and @u/én on the remainder. 


A solution to the problem is then a function и which is continuous on DUC and 
satisfies the equation and the boundary conditions. 


SAQ 5 
Show that the membrane equation in Cartesian coordinates. 
a2 a2 
Cu Cu 
=з + тз = —Р(х.у), 
ах? ду? 


transforms into 


= —G(r, 0) 


in polar coordinates (r.0) where x = rcos 0, у= rsinÜ and G(r, 0) = F(r cos 0. 
r sin 0). 


If G(r, 0) = r inside the circle C:x? + y? = a? and u = 0 on C, find a solution which 
is bounded and has radial symmetry (i.e. depends on r alone). 


(Solution on p. 34.) 
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SAQ 6 


Show that 


2 


(r0) = Zu - Jin 0, 
where (r, 0) are polar coordinates and a, b are constants, satisfies Laplace's equation 
V^v = Oforallr > 0. (ніхт: Use the form of the Laplace operator obtained in SAQ 5.) 


Sketch the curves V(r, 0) = с for various values of the constant c in the domain г > a. 
What do these curves look like for large r? 


In the context of fluid dynamics, is called a stream function and the curves along 
which the value of is constant are the streamlines for incompressible, non-viscous 
flow. The fluid velocity is the vector field represented in rectangular Cartesian 
coordinates by 


(x, y) — AE oet 
` ду дх 


Can you guess what the form for-¥ in this problem represents physically? 


(Solution on p. 35.) 


3.1.2  Laplace's Equation and Boundary Data 
READ W: page 50, line —6 to page 52, the end of Section 10. 


General Comment 


In Hadamard's example small changes in boundary data lead to large changes in the 
interior solution since 


lim (7, y) 
nm 


does not exist for any fixed y # 0. Since lim f(x) = 0 for large n, and u(x, y) = O isa 
solution of the problem in which f is the zero function, we might expect the limit of 
ux. y) for x € (0, 1), y > Oto be zero too, which is clearly not the case. This is what we 
mean by saying that the solution is not continuous with respect to its data. 


A boundary or initial value problem for a partial differential equation is said to be 
properly posed if and only if its solution exists, is unique and depends continuously on 
the assigned data. Generally speaking, a problem which is not properly posed is 
unlikely to represent a physical problem. 
SAQ 7 
Show that 

ет “"соѕһ(4п + 1)лу | neZ* 
is unbounded for any fixed у # 0. 
(Solution on p. 36.) 


SAQ 8 
W: page 52, Exercise 3. 
(Solution on p. 36.) 
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3.1.3 Green's Theorem 


Before you read the next passage from И” we need to remind you of, or introduce you 
to surface integrals and line integrals. The former are described in Unit MST 282 6, 
Rigid Bodies, Section 6.2.3 and the latter in Unit MST 2827.Section 7.1.1 :alternatively, 


you could refer to Spiegel, Vector Analysis, Chapter 5. The essentials are summarized 
below. 


The surface integral i f dA of f over the surface D represents the limit of a sum 
D 


formed by adding together the products fpAA where f, is the value of fat a point P in 
the small element AA of the surface D. When the surface lies in a plane the integral 
may be evaluated by introducing coordinates so that it can be converted into an 
equivalent double integral. In a rectangular Cartesian frame of coordinates Ox). we 
have AA = AxAy and the integral is evaluated as 


. 
[fes nas ay: 
JJD 


In polar coordinates (r, 0), we have AA = rArA0 and the integral becomes 


r 


j f(r cos 0, r sin 0)rdr 40. 
Jp 


(In a proper treatment we should check that these two integrals, are, in fact, equal.) 
To interpret the double integral we must imagine the inner integral [res y)dx, say) 


to be enclosed in brackets. The limits of integration must be suitably determined and 
it is not always a simple matter to obtain them for a given domain D. 


The order of integration, that is, whether to evaluate 


| JG y)dx dy 
D 
J , Oy dx 


does not normally matter, but a judicious choice will simplify the limits of integration. 


or 


Example 


Evaluate ff JdA where f (x, у) = xy and D is the area shaded in the diagram. 
D = 


КВА 


Solution 


In this case 


1 axti 
ff xydA = f if sy dha 
D “0 9 
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(The braces are not normally used; we have incorporated them here to remove any 
doubt as to which integral sign refers to which integration.) 


1 2x41 1 y axel 
| li xy z dx = | |5] dx 
о (Јо o L2.o 


[ xQx + Du 
= a dx 
о 


(It is not essential to integrate over у first— we could integrate over x first. In this case 
it would be more complicated to do so.) 


SAQ 9 


Evaluate 
Í divy dA 
D 


where D is the disc 

пох, у):х® +<) 
and 

у= x/l = xS y/o xS. 
(Solution on p. 37.) 


The line integral [ f ds represents the limit of a sum formed by adding together the 


products of fp As Where fp is the value of f at a point P in the small element of 
length As of the curve C. We use the notation 


f. 
to indicate that the curve C is closed, and unless otherwise stated, that C is described 
counter-clockwise. If the curve is given by 


$i (х, у) тє [to. 1i) 


then the arc length is given by 


*f[dxP [dypl 
ГИТ 


so that 
ds _ dxY | dy? Д 
dv Ца ат 
For example, if is given iri terms of x we obtain the identity 
TN gui 
dx dx 


The line integral f / ds then becomes the definite integral 


IL + (2) | dx 


between appropriate limits. 
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Example 

Evaluate } 
Cc 

and n is the unit outward normal to the circle at the point (x, v). 


v- nds where C is the circle x? + y? = 1 and v = x/ = х?) + wu = х2), 


у= б х3) 


“+ 


у= Vl - xh 


Solution 
The unit outward normal n is given by 
EH EIS 
егай (х2 + y?) 
(see Section 3.1.1). Thus 


xi + yj 


= xi + у] 
since x? + y? = 1 on C. Hence, 
үп = xu = х?) + ya = х?) = Ju = x’) 
on C. Let C, be the upper part of C, on which 
y= Jl - x) 
and C, the lower part, on which 
y= -yU - x» 
On C, we have 
dy x 
dx — JU- x?) 
so that 


ds 
dx 


Ш 
eon 
+ 
—— 
z.|.R- 
21 

о 
ME 


[ 
+ 
! N 


TJU) 
taking the negative square root since we are traversing C counter-clockwise, and so 
ds/dx is negative on C,. Thus 


ae! af 1 
= = x) 57 dx = c= 2. 
NL Ju = x? Jü Ta% f & 
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Similarly, on C; we have 


ds — EX we 
dx y= x?) 
Thus 
" H 7 à 1 1 
Zr рб = | dx = 2. 
k es VI = x?) -1 
Finally ` 


| аве | vndis | v:nds = 4. 
Jc e с: 


It is not by chance that this numerical result is the same as that for SAQ 9; it is а 
particular example of the Divergence Theorem which you will meet shortly. 


‚ READ W: Section 11. page 52 to page 53, two lines below Equation (11.2). 


Notes 


(i) W: page 52, Equation (11.1) . 
Here, div, grad and V? аге to be interpreted as operators on functions of two 
variables. An analogous result holds when the operators act on functions of 
three variables. 


(i) W: page 53, lines 5 to 8 
A proof of the Divergence Theorem is given following this reading passage. 


The name Stokes Theorem is usually reserved for a particular generalization 
of this result to higher dimensions. 
DIVERGENCE THEOREM (IN TWO DIMENSIONS) 


Let D be a convex domain in the plane bounded by the closed curve C, and let 
v = (vy, vy) be a vector field, where (x, у) — v, and (x, y) +> v, are continuously 
differentiable in D and continuous in D u C. Then 


I divvdA — { vends, 
D JC 


where n is the outward unit normal on C. (A domain in R" is said to be convex if no 
straight line intersects its boundary in more than two points.) 
Proof 


Assume that C is bounded by the lines x = a and x = b (a « b), as shown. Let 
y = f(x) on the upper part of C, and y = g(x) on the lower part. 


y ^ 
C y = f(x) 
В 1 
| 4 
| iy = g(x) 
| i 
0 а” = D 
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Then 


Ш 


rr ap b ( ftx 
[| Daa [| i Fea ax 
p Cy TE 


b 
= (вло tea us) dx 
a 


І 


b a 
-f 0-х -Í vl, sey 
a b 


= { vydx. 
С 


А similar argument іп which the integration with respect to x is performed first gives 


дь 
2 ФА = pde 
Í 5 Ox A } v,dy, 


but note that a different partition of the curve C will be required. Addition gives 


ff divvdA = $ v,dy — vydx. 
D c è 


To show that 


j vends = f куйу — v,dx, 
c 


we note from the figure that Ax = —Ascos0, Ay = —Assin 0. (Note that C is 
traversed counter-clock wise.) 


The unit outward normal is 
n= —sin 01 + соѕ0ј 


dy, dx, 
- ds ds 


^ 


so that 


j v,dy — v,dx -4 v:nds 
e c 


If a domain is not convex, but can be divided into a finite number of convex sub- 
domains, then the theorem holds for each sub-domain, and for the whole domain by 
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addition, since it turns out that the line integrals along internal boundaries cancel in 
pairs. (In the figure the contributions along AB and BA balance each other.) 


A B 


Example 


Find the area enclosed by the ellipse 


using a line integral. 


Solution 


The area A enclosed by the éllipse is defined by 


a= ff aa 


(You may verify that this definition of area is equivalent to that in Unit M100 9, 
Integration 1.) Now put v = (x, у) and apply the Divergence Theorem. Since 


divv = 2, 


we obtain 


1 
di 5j (xdy — ydx). 
2с 
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We can represent the ellipse by 
t (a cos г, b sin t) t e [0, 2л], 


the description being counter-clockwise as t increases. (Check that this satisfies the 
equation for the ellipse.) Thus, by the rule for substitution, 


(a cost. b cost + bsint.asin t)dt 


Ш 
$—3 
3 


о 


io» |, dt = nab. 


The area obtained agrees (unsurprisingly) with the result in Unit M100 13, Integra- 
tion 11, Section 13.2.5. 


SAQ 10 


Evaluate f v: nds where v = (—5x, —2y?) and C is the triangle whose vertices are 
e 


(0. 0). (1, 0). (0, 1). Verify that it equals Í div у 4А where D is the interior of the 
D 
triangle. 


(Solution on p. 37.) 


SAQ 11 
Establish the following identities: 
(a) div(a b) = adivb + grada: b: 


Ay: 
I иу? dA -} u £ ds - | grad u: grad v dA; 
D c en D 


5 др à 
| (uV3v — vV?udA = f | и = -v 5] ds; 
D c\ On дп 


where а, и and г are scalar fields and b is a vector field, and C is the boundary of the 
convex domain D. 


(Solution on p. 38.) 


PDE 3.1.4 
3.1.4 Uniqueness 
READ W: page 53, line — 14 to page 54 (the end of Section 11). 
Notes 


(i) W: page 53, lines — 7 and —6 
u, and и; both satisfy Equations (10.1) on W: page 48; that is 


V7u, = —F inD, 


щи =f on C, 
for i = 1,2. By linearity, the difference v = u, — t; satisfies 
Viv20 inD, 
v—0 onC. 


(ii) W: page 54, line 1 
This is a mixed boundary value problem in which values of и and @u/én are 
given over different parts of the boundary. 


SAQ 12 
Show that the boundary value problem 
V?u=0 inD, 
P (x 9) + h(x, Hux, y) = g(x,y) and (х,у) > 0 onC, 


where C is a piecewise continuously differentiable curve bounding the domain D, 
has at most one solution which is continuous on D u C. 


(Solution on p. 39.) 
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3.2 EXTREMUM PRINCIPLES 


3.2.0 Introduction 


It is often impossible to find analytical forms for solutions to Poisson's (or Laplace's) 
equation even when the data are specified on a boundary curve with relatively simple 


geometry. Numerical techniques are required and these are discussed in detail later 
in this course. 


In this section we develop some general inequalities which provide bounds on the 


value of the solution и in the domain D on whose boundary data are given, and on the 
integral 


ff и аА. 
р 


We have set the discussion in the context of fluid mechanics, with particular reference 
to steady flow of a viscous fluid in a pipe. although the techniques can also be used 
with success in analogous problems such as the twisting of elastic bars. Generally 
speaking the method will not give exact solutions, or even good local estimates of 
solutions, but it remains an approach which often provides results when most other 
forms of non-numerical analysis fail. 


3.2.1 The Maximum Principle 


READ W: Section 12, page 55 to page 56, line 12. 


Notes 


(i) W: page 55, line 2 
Although the Law of Conservation of Energy was used to prove the uniqueness 
property, we derived it mathematically ( W: page 53) from the problem. and not 
from a consideration of the physics. 


(ii) W: page 55, line 8 
uattains its maximum on D vu C (denoted in W by D + C), asaresult ofa theorem 
in analysis which says that if a function, with domain a closed, bounded subset of 
К" and codomain R (the reals), is continüous, then its values are bounded and it 
attains its bounds. . 


(ili) W: page 55, lines 9 to 14 
Clearly, since u has a maximum at (xo, yo) € D, the function x > ul,.,, has a 
maximum at xg, and the function y +> u|, -,, has a maximum at yg. Now, we 
know that at a point where a function of a single variable has a stationary value 
its derivative is zero, and for a maximum the second derivative is non-positive 
(Unit M100 15, Differentiation I1). Hence, 


ёи eu 
= ыйы < 
= =0, a <0, 
OX oso) X^ [xo.yo) 
and 
" 
си 
РТН = 0, 
CY [Go yo 
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maximum 
value 


if 


On C, however, it is possible to have a maximum which is not a stationary value. 


(iv) W: page 55, line 18 
Note that the foregoing proof must fail if F takes the value zero at some point 
in D since 


2 
V7Ulsxo.yo) © 0 


апа 
Е(хо, уо) <0 


аге not necessarily contradictory. Thus a more elaborate proof is required for 
the case F(x, y) € 0. 


General Comment 


Although we have not asked you to study the proof of continuity with respect to the 
data for the Dirichlet Problem, you should be familiar with the result. 


SAQ 13 


Two harmonic functions (x, y) э u, and (x, y) +> и, satisfy u; < t; on a piecewise 
differentiable closed curve C. Show that и, < u, throughout the domain D bounded 
by C. 


(Solution on p. 39.) 


3.2.2 Flow in a Pipe 


Consider the steady flow of an incompressible viscous fluid along a straight pipe of uni- 
form cross-section. By steady flow we mean that the speed of the flow w is independent 
of time at each point. Incompressibility implies that w must also be independent of z, 
the distance along the pipe (as shown in the figure), for otherwise a given mass of 
liquid would be compressed irito a smaller volume or expanded into a larger volume. 
If D is a section of the fluid in the pipe then it can be shown that w satisfies 


Ур = —k inD. (1) 


Here k = p/pvisa positive constant; p is the (constant and) uniform pressure gradient 
along the pipe, is the density of the fluid and v is its kinematic viscosity, a measure 
of the internal friction of the fluid. A viscous fluid adheres to the boundary (think of 
oil flowing down a metal plate, or why the engine oil level can be measured with a 
dipstick) so that 


w=0 onC. (2) 


PDE322 


This problem is equivalent to that of a membrane, attached to a plane frame, with 
one side subjected to a uniform pressure (W: page 48). 


By the maximum principle of Section 3.2.1, w must attain its minimum on C, since 
Vv7(—w) = k > 0. Thus w > 0 in D; in other words, local flow against the pressure 
gradient is impossible. 


The total flux, defined by 


Q= Í, wdA. (3) 


is of physical interest since it measures the rate of flow through the pipe. Since 
1 = —(V^w)/k, we can express Equation (3) as 


Q- ~ iff vm dA, 


1 ew 1 A 
== i "3n ds + n Jf ioa wl? dA, 


using Green's Theorem (W: page 53). However, the first integral on the right vanishes 
by Equation (2). Thus 


Q= All lgrad wl? dA. 
к 


SAQ 14 

Verify that the fluid velocity w in a circular pipe of radius a is given by 
w = 4k(a? — x? — y?) 

where k is as defined in the text. Find the flux through the pipe. 


(Solution on p. 39.) 


SAQ 15 


Find an expression for the speed of flow in a pipe with an elliptic section given by 


2+2=1, 
at BR 


when К is as given in the text. 


HINT: Guess the solution by generalizing the result for a pipe with a circular section, 
and then verify. 


PDE 3.2.2/3.2.3 


Show that the flux is given by 


3.2.3 Extremum Principles for Viscous Flow 


In order to derive bounds on the flux for the viscous flow problem, we first obtain a 
result about vector fields on a domain D in R?, that is, continuous functions D —> G?, 
where G? is the set of geometric vectors in the plane. 


We have seen in Section 3.1.1 that the set of vector fields on D forms a vector space 
with respect to pointwise addition and scalar multiplication of vector fields. 


For two vector fields 
a:D—>G*, рэб. А 


let a - b: D — R denote the scalar field whose value at each point of D is given by 
the usual inner product for geometric vectors. Thus a -b is the mapping 


(х, у) э а(х. у)" Ыбх, y) (х, уер. 


Now. define the binary operation Ооп the space of vector fields by the surface integral 
(expressed in function notation) 


anb = ff a». 
D 


It is not difficult to see that for any two vector fields a, b on D the operation O yields a 
real number and satisfies the axioms for an inner product (you need not verify this 
statement). With respect to this inner product, Schwarz’s inequality (see Appendix) 


ТЕТЕП 


We shall employ this result to рш bounds on the flux О, which, as we have seen in the 
previous section, can be expressed as 


id, jgrad w|? aA (4) 


for viscous flow in a uniform pipe. 


FIRST EXTREMUM PRINCIPLE 


Let w satisfy V^w = —k in a domain D and w = 0 on its boundary C. Then the total 


flux 
Q= fi wdA 
D 


is bounded below by 


Ф 
k І w* aa i lgrad w*|? dA, 
D D 


where w* is any non-zero scalar field on D ùu C with the property that w* = 0 on C. 


PDE 3.2.3 
Proof 


If we set u = w*, v = w in the identity proved in part (b) of SAQ 11, we obtain 


ff w*V^wdA = { we == ds - | grad w* - grad w dA. 
D D 


c n 


Since w* = 0 on C and V^w = —k in D, this reduces to 


k Ij w*dA = ff grad w* - grad w dA. 
D D 


We now apply Schwarz's inequality with a = grad w* and b = grad w, to obtain 


e(ff w* aa) І < [| lgrad w|? dA || [ргай w*|? dA 
D D D 


x kQ ff lgrad w*|? dA using expression (4) for Q. 
D t 


The result follows since k > 0, and so 


Q2k (T. w* da) lif, lgrad w*|? dA. 


SECOND EXTREMUM PRINCIPLE 


Let w satisfy V^w = —k in a domain D and w = 0 on its boundary C. Then the total 


flux 
Q= ff wdA 
D 


is bounded above by 


1 
t S 
il, v** dA, 


where v* is any vector field which has the property that div v* = —k in D. 


Proof 
Since div(wv*) = w div v* + v* : grad w at all points of D, by part (a) of SAQ 11, 


ff v*- grad w dA = I div(wv*) dA - | w div v* dA 
D D D 
= { wv*- nds — ff w div v* dA, 
Jc D 


using the Divergence Theorem (Section 3.1.3). Since w = Ооп Cand divv* = —kinD, 
this reduces to 


] v* - grad wdA = КО. 
D 


by the definition of Q. Applying Schwarz's inequality with a — v* and b — grad w, 
we obtain ] 


EQ < ff |grad w|? dA ff v*? dA 
D D 
< ко ff v dA 
D 


using expression (4) for Q. 


The result follows since k > 0. 


M 
o 
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Fields w* and v* which satisfy the conditions above are said to be admissible functions 
for the First and Second Extremum Principle, respectively. We have established that, 
for admissible functions w* and v*, 


Le oy 
nmm [гай w*[ dA «iJ. 


Example 


Find bounds for the flux in a pipe of square section, (x, y) e[- a, a] x [—a. a]. 
i 


Solution 


To obtain the upper bound we choose v* = (—4$kx, —4ky) which is admissible since 
it satisfies div v* = —k. (This could be used for a pipe of any cross-section.) Thus, by 
the Second Extremum Principle, 


[os nf f (x? + y?) dx dy = Ska*. 


For the lower bound we require a function which vanishes on C; we choose 
w* = (a? — x*)(a? — y?) 
Thus, by the First Extremum Principle, 


di, ie (a? — х?)(а? -dxd | 


шя fia [x a? — y? + уа? — x*y)dx dy 


= ёка“. 


(Verify the integration if you have the time.) 
We conclude that 
0.555ka* < Q < 0.667ka*. 


The correct value for Q, to three significant figures is 0.562ka*. (This result will be 
obtained from the eigenfunction solution to the problem in Unit 6, Fourier Series.) 


The construction of the functions w* and v* requires some practice to gain a feel for 
the most useful choices. It might help to give some idea of the thinking behind the 
choices in the example. In general, keep your choice of functions as simple as possible. 
For w* you must select a function which vanishes on the boundary. Thus, for the 
square pipe we chose 


w* = (a? — x?)(a? — у?), 


which vanishesonx = +aandy = ta. If we think of the analogy with the membrane 
problem, we can visualize the shape assumed by a membrane attached to a square 
frame when it is subjected to uniform pressure on one side. We would expect the 
maximum displacement to occur at the centre; this also occurs for w*. 


For equality to obtain in the upper bound for Q we should have v* = grad w. Now, 
in the case of the membrane analogue we would expect, by considerations of sym- 
metry, that ди/дх = 0 on x = 0 and éw/dy = 0 on y = 0, so that ož = 0on x = 0 
and v = 0 on y = 0. The choice of v* = (vf, v) in the example is the simplest one 
which has this property in addition to satisfying 


divv* = —k. 


The following SAQ is lengthy, particularly part (ii). It is incorporated because it 
enables you to see how accurately the extremum may be determined. 
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SAQ 16 
(i Show that 
y* = [—Жх + p(x? — 3xy?), —dky — n(3x?y — у?)] 
is an admissible function in the Second Extremum Principle (i.e. divv* = —k) 
for any value of the constant д. 
(i) Show that the corresponding upper bound for О in the square channel of the 
example is 
a*[2k? — l6kua? | 96д?а* : 
К\з is * 35 


(iii) Determine the value of и which minimizes the upper bound and find the least 
upper bound. 


(Solution on p. 41.) 


Here is a similar problem you can try if you have the time. 


SAQ 17 


Show that v* = (ux, (—k — j)y) is an admissible function for the Second Extremum 
Principle for any и. Show that the best и gives an upper bound for the ellipse 


En 
do 


equal to the exact result obtained in SAQ 15. 


(Solution on p. 41.) 
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3.3 THE DIFFUSION EQUATION 


3.3.0 Introduction 

The partial differential equation which provides a good model for time-dependent 
heat conduction is also applicable to diffusion, e.g. 

(a) diffusion of mixtures of liquids or gases, 

(b) diffusion of vorticity in viscous fluids (vorticity is the local spin of the fluid). 


The one-dimensional form of the equation for heat conduction was derived in Unit 
M201 32, The Heat Conduction Equation. In this section we: shall indicate how this 
can be extended to three dimensions. Before we do that however we need to develop 
the three-dimensional form of the Divergence Theorem. This in its turn requires some 
familiarity with non-planar surface integrals and volume integrals which we therefore 
introduce. 


3.31 The Divergence Theorem 
In Section 3.1.3 we introduced the concept of an integral over a surface contained in 
the plane. For surfaces in three dimensions we adopt the following approach. 


A surface in R? is defined as the image under a one-to-one map whose domain is a 
connected subset of R? and whose codomain is ЁЗ. It is often convenient to represent 
the points of the surface by position vectors in КЭ. (As for curves, we usually require 
the surface to be piecewise continuously differentiable, but a precise formulation of 
this condition for surfaces would take us too far afield.) 


Thus, 
S = {r:W(u, v) = r, (u v)e M) 
is a surface for a suitable M c R?, and г = (x, у, z) e R?. For example, 
r = (a cos u sin v, a sin u sin v, a cos v) u € (0, 27), v e [0, 47] 
defines a hemisphere of radius a, and is the image of a rectangle in the ur-plane. 


v 


If the functions u —» x, и» y, wz, v > x, о э у and р> 2 are con- 
tinuously differentiable at the point P whose position vector is r — (x, y, z), then we 
define the vectors 

er Ox ду дт 

ди ди ди ди 


апа 


or Ox ду Oz| 
др 
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The arrows belonging to these geometric vectors which begin at P are tangents to the 
curves through P along which г and и arc constant, respectively. Together they span 
the tangent plane at P (i.e. the plane consisting of arrows at P which are tangent to the 
surface at P). The vector 


п (ey 00] 
BU" ejl 


is then a unit vector normal to the surface. 


т а 


X ста 
ёи ёг 


We shall now look at a neighbourhood of P in more detail. We shall use the function yr 
to define a coordinate system on $: thus we may say that the point P at which 
= (и. v) has coordinates (u, v). 


(u + Ane) 


0 


The surface is covered by a mesh of the two families of curves determined by и = con- 
stant and v = constant. Let P, be the point y(u + Au, v) and Р, the point (и. v + Ac). 
The vectors representing the chords PP, and PP, are given, in the first-order Taylor 
approximation, by 


PP, = W(u + Auc) — Wu, v) = T Au, 
— ди 
РР, = y(u, v + Av) —1Д(и,р) = ль. 
=$ ёр 
Thus, to the same approximation, the area of PP,P,P, is (see Appendix) 
AS =|PP, x РР, = [< x Claude. 
—> ^ —$ Ou дг 


This represents an increment of area of the surface. To integrate a scalar field over 
the surface S we define the surface integral- 


ff pds 
s 


Jm 2 #48 


[А 


Thus any surface integral reduces to an integral over a surface їп the plane which can 
be evaluated as a double integral over Cartesian coordinates. (Where convenient, the 
double integral in the plane may be transformed to other coordinates, e.g. polar 
coordinates.) Although we have used S to denote both the surface and an element of 
area, this should cause no confusion. 


or 
х — 
Ov 


er 


du dv. 


The surface area of S is defined as [| dS: as above this equals 
5 


UNE 


= x C [du dv. 


ди ёр 
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The volume integral fif dV represents the limit of a sum formed by adding 
D 


together the products $5A V where фр is the value of ġ at a point P in the small element 
AV of the domain D in R? (see Unit MST 282 6, Rigid Bodies, Section 6.2.2). We 
evaluate such an integral as the triple integral over rectangular Cartesian coordinates 


Iii фах dy dz, 


between appropriate limits. 
We now have sufficient machinery to state the Divergence Theorem in R?. 
DIVERGENCE THEOREM (IN THREE DIMENSIONS) 


Let D bea convex domain in R? bounded by the closed surface S. Let the coordinates 
of the vector field v have continuous first partial derivatives in D. and let n be the unit 
outward normal to S. Then 


fff iva = ff nonas. 


Proof 


E 


Divide S into two surfaces S, and S, as shown and assume that they can be repre- 
sented by z = f(x, y) and z = f(x, y). Let v = (v, , v2, v3) and consider 


zcfiny) 
eo Leelee 
м =» 92 


where M is the projection of D on the xp-plane, i.e. М = z(D) where л :(x, y, 2) — (x, у). 


Thus 
д 
ff SV = Jf. [os] Z5 dx dy 


The surface S, is given by 
r-(s»fioy) (x y)eM, 


and so 


ôr _ (1.04), дг _ [o T 
Ox 
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and 


The outward unit normal on S, is 


Ж т cr| Даг er 
= % =} | 
Ox  Oy][|C 


X 
Ox ду 
The minus sign arises because the z-coordinate of n must be negative. By definition, 


I vk: ndS = [], v3k -n|— x x 
Sı 
- -ff vk: (= « on - 


dx dy 
ôx ду 


Now, 


and therefore 


I u3k-ndS = -ff 03|. - p s dX dy. 
Sı M 


A parallel argument gives 


| vak ‘ndS = lj 03|. - pae dx dy, 
3; M 


and thus 


Wa uw ff vak : n ds. 


Similarly we may show that 


iff Sav = ff v,i n dS, 
p 0х 5 
if. = dV = I vaj’ n 45. 


Ov, | Ov, дь, 
div v = — + = + = 
WY = x Oy 02 


Since 


and 
v = pb + 05) + vak 


we obtain 


[аа = || пав 


as required. 
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Example 

Verify the Divergence Theorem when v = (xyz, .—) and D is the domain 
bounded by the planes x = 0. x = 1.2 = 0, у = 0. y +2 = 1. 

Solution 


We first observe that 


divy = X (xyz 9+5 Tears (— yz?) = 3x? yz. 
сх ez 


The domain D is shown in the figure, and the vertices are labelled. To integrate over 
the whole volume we integrate first with respect to z, between 0 and 1 — y, then with 
respect to y between 0 and 1, and finally with respect to x between 0 and 1. Thus 


ЇЇ, Е i ARIA “xyz аја 


1 1 
= з}, [ 4x71 — yy*dy dx 


-if ухх = зц. 


We now have to evaluate [ v:nd$ over the surface of D. First note that the unit 


outward normals n on the respective faces (which are denoted by P, Q, R, S, T) are 
as follows. 


у= 0:пр = –ј 

z = 0:по= – К 

х= 0:пр= 7i 

х= 1:0; = і 
y+z= l:nr=(j + k)//2 


We thus obtain 


[f rome = [f v reas = ff v nas =o. | 
IE +0545 = ACC yz dlas оа 


Jf ‘ny dS = zh [ (y?z - уг?].-1-уу/24х dy 
E: f [ DPU = y) = y — у)914х4у = 0. 
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Thus, we see that 


fij. div v dV = 30 = [| y-nds. 
surface of D 


SAQ 18 


Verify the Divergence Theorem when v = (xy, yz, zx) and D is the cube bounded by 
x=OQx=Ly=O0y=l2=02=1. 


(Solution on p. 42.) 


3.3.2 Derivation of the Heat Equation 


Suppose the material in which heat is flowing occupies the domain D c R? and let 
u(x, y, 2, t) be the temperature at the point (x, y, z)e D at time t 2 0. Fourier's heat 
transfer law is 


J = —K gradu, (1) 


where J is called the heat current vector and represents the rate and direction at which 
heat is flowing at a point, and K is the thermal conductivity of the material which is 
assumed (for our purposes) to be constant in time and space. Fourier's Law states that 
heat flows in the direction of greatest change in temperature from points of higher 
temperature to points of lower temperature (this accounts for the minus sign). 


Consider now any subdomain Dy of D. The heat content H(t) of Dg is defined to be 


H(t) = ср ЇЇ, u(x, у, z, t)dx dy dz, (2) 


where p is the density and c the specific heat of the material; both are assumed to be 
constant in space and time. (The specific heat is the heat required to raise the 
temperature of unit mass of the material by one unit.) The heat equation follows 
from the Principle of Conservation of Heat (a form of the Law of Conservation 
of Energy) which states that the rate of heat accumulating in Dg equals the rate of 
heat entering through the boundary; in mathematical terms 


—- а J-nds, 
dt Co 


where the surface Со is the boundary of Do and n is the unit outward normal to Cy. 
Thus, from Equations (1) and (2), and differentiating under the integral sign, 


cp [if = dx dy dz = K | (grad и) > п dS 
Do бї Co 


=K ff V?u dx dy dz 
Do 


by the Divergence Theorem. Finally 


I { o% - KV? uh dx dyd = 0; 
l.. 


since this is true for any domain Dy we infer that the integrand vanishes, giving the 
heat or diffusion equation 
ĉu 
c = КУ?и, 
д 


where К = К/ср. 
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If the temperature is steady, then и is independent of t and must satisfy Laplace’s 
equation 


VÀu = 0. 


There are three main types of spatial boundary condition which may be satisfied by 
the temperature: 


(i) w is a prescribed function of t at each point on Co; 


(ii) n-gradu is prescribed on C, (the boundary is said to be insulated at those points 
at which n : grad u = 0); 


(iii) n- grad u + (и — ug)h = O where h and ug are constants (this represents radiation 
from the surface into a medium of fixed temperature). : 


3.3.3 Uniqueness and Maximum Principle 


READ W: page 59, line 1 to page 61, line 5 


Notes 


(i) W: page 59, lines — 12 to —7 
This represents a statement of the uniqueness theorem for u = и, — u, where 
u, and wu; are any two solutions to the heat equation in D such that u, = u, on C 
for t e [0, t] and in D for t = 0. It is then shown that u = 0 іп D for 0 x t < t. 


(i) W: page 59, line —4 
Green's Theorem in three dimensions is a special case of the three-dimensional 
Divergence Theorem (Section 3.3.1) with v = u grad u. 


SAQ 19 


Show formally that, for any function f: R —+ А which is continuous and whose 
values are bounded, 


© 
u(x,t) = xu] g^ 6-3? HM Fe) dE 
-® 


2(ktn)i 
satisfies the one-dimensional heat equation 
Qu д?и 
= К = 0 
ді ox? 


(You should assume that you may differentiate under the integral sign.) 
By substituting č — x = 2n(rk)* in the integral show that 

lim u(x, t) = f(x), 

n 


under the assumption that, for the integrand concerned, 


" " 
lim | - | lim 
1-0 Ju 2179 


© 
mnr: | er dy = n. 
-0 


Describe a heat conduction problem to which the function given above is a solution. 
(Solution on p. 43.) 

SAQ 20 

W: page 61, Exercise 2. 

(Solution on p. 44.) 


32 


PDE 3.4 


34 SUMMARY 


We have covered the following main points in this text and the associated reading 
passages: 


(a) uniqueness of solutions of Poissor's equation and the heat equation with appro- 
priate boundary conditions; 


(b) the maximum principle for the same two equations; 
(c) the notion of a properly posed problem. 


The application of extremum principles to flow in a pipe was introduced and we 
showed how bounds could be obtained for the flux through the pipe. 


The following results were introduced: 


Stokes' Theorem in the plane, 
the Divergence Theorem in two and three dimensions, and 
Green's Theorem in two and three dimensions. 


In the course of discussing these theorems we defined curves and surfaces in R? and R? 
and developed techniques for evaluating line integrals, surface integrals and volume 
integrals. 
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35 SOLUTIONS TO SELF-ASSESSMENT QUESTIONS 


Solution to SAQ 1 


(iii) 0 

(iv) 3j — 3k 

(у) 6 

(vi) 6i — 3j — 3k 

Solution to SAQ 2 

(а) 3(x? + y? + (xi + yj + zk) 
(b) yzi + zxj + хук 

(c) 2xi + 2yzj + (у? + 322) 
Solution to SAQ 3 

div F = 2y + 6xz. 


Solution to SAQ 4 


A unit vector in the direction given is 


At (1, 2, 1) 
grad $ = 2i +j + 2k. 
Therefore, at (1, 2, 1) 
дф 1 1 
= = е: ргайф = —= (2-1 - 2) = –- =: 
ĉe J^ "E 
Solution to SAQ 5 
By the chain rule for functions of several variables, we have 
ди ди Ox | Ou Oy 


or Ox Gr дуд 


and 
Ôu Gu Ox ĉu? 
A epee peram 
д0 х д0 у 80 
Since 
x=rcos0, y=r sin 0, 
we have 
Ox E , ] 
X = cosh =~, 2 = sinf = 2. 
дг “д " 
and 
Ox 
26^ —rsinü = —y zz = rcosô = х. 
Thus, 
ĉu ĉu ди 


à хах 72у 
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and 
ди ёи ди 
E ы oce betel 
cü 7 дх ду 
Непсе, 
êf ди ё " 6 ди Е ди 
r—|r =(x-= y x y—]> 
ér\ дг ex  "OQyJV ex ду, 
ie. 
E д?и HL „д?и ди eu au 20u n ё 
. r—=x x yx | , у= 
ĉr? r ex? ex дудх дхду ey! "ày 
and 
Uu д " д ди R Qu 
={-y x у х 
80? оу (7а ар 
eu ди eu д?и ди eu 
=y 1T-xXi-xXyz-—gXi—y х?-— 
"ox ex 7 дудх дхду ду ду? 
Therefore, 
д2 a2 2 a2 
,0?u ĉu Qu A a[O д?и 
PP tra + Spa? D 
or? er 00? ( ) Ox! ду? 


and the result follows, since r? = x? + y?. 


Let w:r — и; then, with G(r, 0) = r and assuming that u has radial symmetry the 


equation becomes the ordinary differential equation 
м (т) + two) =-r 

or 
КЕ? =-r, 


Integration gives 
rwo) = =} + А, 
where A is an arbitrary constant. A second integration gives 


w(r) = —i? + Alnr + B, 


where B is a further arbitrary constant. Since In r is unbounded as r ~ 0+ (i.e. as r 


tends to zero through positive values), we put 4 — 0, whilst B is given by 
0 = ҙа? + B. 
Thus the solution is 


и = (аз — r°). 


Solution to SAQ 6 
The derivatives of 'F at (r, 0) are given by : 


Mir 0) = of + 5 sin 0, 


er 
eg) —2ba^ . 
a (0) = 3 sin 0, 


and 
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thus 
ep 109+ 1 ay 
2 tra RW 
The figure shows some typical curves along which the value of V is constant near the 
circle r — a. 


= —_—_ 
=_= 
= = 
с С = 
= = = 
== = = 


For large г, 
Wir, 0) = br sin 0 = by. 
Thus the curves approach the straight lines y = constant. 


In the context of fluid dynamics, the velocity field for large r is approximately (b, 0), 
a uniform stream. The curves in the figure represent streamlines; note that 
Wr, 0) = V(r,z) = V(a,0) = 0 for r > a and 0 € 0 < 2л. 


Our choice of ¥ represents uniform flow past a circular cylinder. 


Solution to SAQ 7 

e^ "" cosh(An + 1)лу = е7 "(e^"* 0n» — go mt m) ; 
For large n and y # 0 this behaves like 

^nm i 


For 


n > —- 
16r? y? 


this has a positive exponent which increases without limit as n becomes large. Hence 
the expression is unbounded for large n. 


Solution to SAQ 8 
div I = div(cE) = — div(c grad u). 
Thus u satisfies 


div(c grad u) = 0, 


e s” n д{ ёи ef ди 0 
24° 2х) + ду\ dy} "ad az] = 


до ди Е до Qu до ди 
дх дх дуду! Oz dz 


ie. 


or 


oV7u + 
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The second derived functions occur in the term ¢V?u which is the product of (0) 
and the Laplace operator acting on и. Thus the equation must be elliptic. 


Solution to SAQ 9 


divv = Že АЛ x?) +20759) 


" 2 — 3x? . 
Ju — x?) 


у= Jf(l - x) 


ny 


у= -JS(l - x) 


va 7x2) 
divvdA = f. Г. = dx 
Jf, NE zs 
l - (17x32) 
(2 — 3x?) [ " 
——À |y dx 
alt=" |. аа 
А ; 
= | 2(2 — 3х?)ах = 4. 
-1 
(Integrating with respect to x first would have been rather tedious!) 


Solution to SAQ 10 


B(0, 1) 


[v A(LO) x” 


1 
Í v:nds = | [257], odx = 0, 
OA о 


since n = —j and у = 0 on OA. 
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0 
I v:nds = | [5x].-ody = 0, 
BO 1 


since n = —i and х = 0 оп BO. 

The line AB has the equation x + у = 1. In this case the line is given parametrically by 
x=l-y ye (0, 1]. 

The outward normal is n = (i + AE and 
е 
dy dy 


= fi. | 


Therefore 


iL 
f yends =f [-5(1 — y) — 2y?]dy 
AB 0 


bp ucl 
1 ix д д 
| divy dA = [ f Е (= 5х) + = (25?) | dy dx 
D o Jo ex ду 
1 1-x 
- | | [—5 – 4y)dy dx 
о Jo 


- f [-5(1 = x) - 2(1 — х)?]4х 


0 


li 


1 
| (—7 + 9x — 2х?)4х 
o 
-— 
2, 


Solution to SAQ 11 


Q(ab,) | ĉlab,) nece А 
ax + ET > where b = bj + b,j 


да ab, да ôb, 


um *3 „жал, 
= adivb + grada b. 


(а) div(ab) = 


(b) The Divergence Theorem gives 


i div(u grad v) dA = f u grad v: n ds 
D С 


др 
= ф и — 05. 
с дп 


Using the identity div(ab) = a div b + grad a: b in the left-hand side: 


[| grad u- grad v dA + [ uVivdA = $ út is: 
D D c On 


Note that if и = v we obtain Green's Theorem (W: page 53) as a special case. 


(c) Interchanging u and v in part (b) we obtain 


I grad v: grad u dA + {I vV^udA -{ v% ds. 
D D c On 


Subtracting this from the result of part (b) we obtain the required result. 
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Solution to SAQ 12 


Let u, and u, be two solutions, which satisfy Laplace's equation and the given 
boundary condition. If v = u, — uz then v satisfies the homogeneous problem 


Vv-0 inD, 


av 
эң w= 0 on C. 


We apply the identity (11.2) on W: page 53 for v, so that 


д 
-$ v OS аз + [f lgrad v|? dx dA = 0. 
c дп D 


Substituting for до/дп on C gives 


{ ho?ds + ff lgrad v? dA = 0. 
c D 


Since h > 0 on C, the integrands are strictly positive unless v takes the value zero 
throughout C and is constant in D. By continuity, v is the zero function on Du C 
and зону = и. 


Solution to SAQ 13 


Let u = u, — uy. It follows that V?u = 0 and u > 0 on C. The minimum value of u 
must occur on the boundary and further 


min u > 0. 
с 
Thus 
uz min и> 0 


in D. Thus u, < из in D. 


Solution to SAQ 14 


If 
w = 4k(a? — x? — у?), 
then 
aw aw 
Ox? %, ay? — ak, 
Thus V?w = —k, so that w satisfies Equation (1). It is easily shown that Equation (2) 


is satisfied, since x? + y? = a? on C. The expression for w gives the velocity since the 
solution is unique. š 


The flux is given by 


epe Io] 


in polar coordinates, where w = 1k(a? — r°). (For our problem, integration is easier 
in polar coordinates than in Cartesian coordinates, since the domain is a circle.) 


Thus the flux is ] 


2n а 
Q= xf \f (ra? — Г, 
o Wo 
= 4kn [ye - | 
o 
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PDE 3.5.0 SAQ 15 


The solution for the circular pipe (SAQ 14) suggests (if you are good at guessing) a 


solution of the form 


which clearly satisfies 


the boundary condition (Equation (2). It also satisfies 


Уу = —k if 
1 1 
-2c|-A + 5] = =. 
Ў Р Ы | 
or 
k db 
РЕЧЕ 
So the speed of flow is given by 
ok ab? i x oy) 
ur a b 


The flux is given by 


bul = x?ja?) 


Q 


ll 


el 
Tl 


2 
he 
| Al 


Ш 


P 


2 
5 - z) o) dx 


5 
I byl 7 x?ja2) 


(1- 


» 


3b? 


(1 x?/a2) 


| 


1—— 


х2 
E dx 
a у= - by T x2/a2) 


Now, if we make the substitution 


х = аѕіп и, 


x= af 
-nz[2 


cos?u . cos u du 


a г"? 


8 |. п/2 


а : ; 
913" + 2sin 2и + 2 sin 4и] 


= (3 + 4 соѕ 2и + cos 4u)du using the hint 


2 
nf2 


= Зла/8. 


mab 


so that 
dx 
— = acos li, 
du 
then 
a x? + 
1- +] d 
Therefore 
4bc Зла 
без" 
_k аЬ? 
“2 e+e 
E пка? 
7 а? +B) 


as required. 
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Solution to SAQ 16 


() divv* = —3k + 3ux? — 3yy? — $k — 3ux? + 3uy? 
= =k. 


and so v* is admissible. 


The upper bound on Q for the square channel is 


(ii) 
di. uae | f {[— kx + ue) — 3xy?)}? +[ iky - щ3х?у— y)P} dx dy 


| ге pe 
=: f f [52 x? – ku(x* — 3x? y?) + p? (xi — 6x* y? 9x? y*) 
+ hk? y? + ku(3x?y? — y*) + u(9x*y? — 6x?y* + y dx dy 


| ге ге 

=i f f [2 (2 + p?) — ku(x* — 6x2 y? + y*) + и(х° + 3x* y? + 3x2 y* + у°]ах dy 
Dm 3 3a5y2 

zi [epe |- nS aee | eae [^ E tay +ay "| 
Г „[а а a$ 24% а gia? at га ае 
| +з uu bas XE dr AE ME 


- a фа +48kna? + 38a) = plu) say. 


(iii) The best и occurs where р'(д) = 0 so that p(x) is a minimum, i.e. where 
pa? = —7k/36. 
Thus the least upper bound is 
ka*[$ — 1$: зв + 38° (55)°1 = тэзКа* = 0.563ka*. 
Compare this upper bound with the correct value given in the text. 


The choice of admissible function in the text corresponds to д = 0. 


Solution to SAQ 17 
v* is admissible since 

divv* = p — k — p = —k. 
The upper bound is given by 


a byti = х20а2) 
iff, v*?dx dy = if if Гах? + (k + ires 
=a (9 —by(1 — х2/а2) 
" by (1 xia) 
= if С + К + й» dx 
ag = by = х27а2) 


) pa xh x2\3 
EE api X Le ap. х r 
х oft =| + 3(k + и)*Ь |: z) |> 


We use the substitution x = a sin u to obtain 


2а b aoe 
ie у= — - Иа? sin?u cos u 4- 3(k + i? b? cos?u]cos u du 
i $ 
-nj2 
2ab pri? 1,242542 1 252 2 
= [qu a sin 2u + z;(k + u)*b*(1 + cos 2u)*]du 
~n/2 
2ab ("^ yg 1 252, 
=T [šu a (1 — cos 4и) + 3; (k + u)*b*(3 +4 cos 2u + cos 4u)]du 
-7/2 
me 


~ urat CE uy! b^]— ри), say 
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The least upper bound occurs where p'(u) = 0, that is, where 
ра? + (К + Б? = 0. 
Thus 
p= —kb?/(a? + b?) 
and the least upper bound is 
2abnk| — a?b* b?a* na^ bk 
8 lc +b gx e] да? +5) 


which is equal to the true value for Q (see SAQ 15). 


Note that if we are interested only in the best value of и we can proceed directly using 
the value of w obtained in SAQ 15. For equality in the upper bound we require 


kb? ka? | 


j* = grad w = | — —£——4 X, — —4——3J 
vB | a? + b? "au 


= (их, (=k — gy) 
when 


и = = Ка? + b?). 


Solution to SAQ 18 


lplel 
I ИИ) (x + y + z)dx dy dz 
D ооо 


ipi 
=f] ($ + y + 2)4у dz 
о Jo 


1 
= [ (5 +44 2)dz 
° 


жаз 1 1-3 
=т=: 


х 


We now evaluate I v: ndS over the six faces. 
s 

For x = 0:n = —i, v = (0, yz, 0), fes = 0. 

For y = 0:0 = =i, v = (0,0,2x), [{v-nds - 0. 


For z = 0:n = —k, v = (xy,0,0), [ча = 0. 
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1 1 

For x = i:n = i, v = (y, ус, 2). fos f Í ydydz = i. 
o 40 
1 H 

For y = n2 j v = (x, 2, zx) fos = f | zdzdx =}. 
ovo 


4 1 
En kv = pao, аа = f Í xdydx = 5. 
о Јо 


ПЕЕ ЕТТЕ fff mar. 
5 р 


Solution to SAQ 19 


For z 


Thus 


Assuming that we may differentiate under the integral sign, we obtain 


ди 1 f | 1 n č- | ет аку (г) dë, 


Ot ^ 2(kn)! 728 ake 

ди 1 © (E — х) а-аа А 

x7 Xa J _„ oe © Sede, 
and 

^u 1 œ T(E — x)? 1 ks 7 " 

sr xe | Saar ало: 
so that 

ди д?и 

a ae 


The substitution 
ë х = 2n(rk)i 


gives 


u(x,t) = wF. ет" f(x + 2n(tk)3) . 2(tk)* dn 


1 fe 
= | ет" (х + 2n(tk))dy. 
Ja | = 


Thus under the assumption that the limit and the infinite integration may be inter- 
changed, 


1 Y = 1 ° =; " At 
lim u(x,t) = AL. e lim f(x + 2n(tk)*)dny 


= fe) E e^" dn = f(x). 


Jalia 


The solution gives the temperature distribution u(x, t) in an infinite bar with a given 
initial temperature distribution u(x, 0) = f(x) for xe К. 
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Solution to SAQ 20 
Extend the domain to [ —1, 1] x [0,1]. as shown in the figure, and let 


u(x, t) = u(— x. t) (x, ))e[-10] х [0,1]. 


1% 


e E 


xy 


0 


Then x +— u(x,t) is an even function and is seen to be twice continuously dif- 
ferentiable over the extended domain, while at the same time 

ди 

(0,0 = 0. 

дх 
In addition и satisfies the differential equation in the extended domain. By the 
maximum principle (in one dimension) the maximum must occur on either t = 0 or 


x = +l. Since u(l, t) = u(—1, 0) we have the required result that the maximum of 
u for (x, t) e(0,7] х [0,1] occurs at t = 0 or x = l. 
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3.6 APPENDIX: VECTORS 


We have collected, in this appendix, some results about vectors which have appeared 
in previous courses and are relevant to this course. This appendix may be used for 
revision or reference. The results are presented in detail in Unit M100 22, Linear 
Algebra I, Unit M201 1, Vector Spaces, Unit M201 16, Euclidean Spaces and Unit 
MST 282 1, Basic Tools. 


A vector space V over a field F of scalars is a set of elements with an internal binary 
operation +:V x V —>» V and an external binary operation of multiplication 
Fx V—+ V satisfying the following axioms for all a,b, ce V and m, ne F: 
1 (a+b)+e=a+(b+c) 
2 at+b=b+a 
3 30eVsuchthatO+a=a VaeV 
4 J-aeV such thata + (—a)=0 
5 m(a + b) = ma + mb 
6 (m+ п)а = та + па 
7 (тп)а = m(na) 
8 la=a 
When F = R (the reals) we talk about a real vector space; when F = C (the complex 
numbers) we have a complex vector space. 
An important notion in a vector space is that of linear independence. A set {V,,.-., Vaj 
of vectors in V is linearly independent if 

«үү, + QV, tb xv = 0 (a; € F) 


holds only for а, = 2; = ::: = a, = 0. Suppose now that {v,,...,¥,,} is a maximal 
linearly independent set in V, i.e. if another vector in V is included in the set we obtain a 
set which is linearly dependent. Then {v,,...,¥,,} is called a basis for V, and every 
other basis for V has precisely m vectors. We say that V has dimension m. It can also 
be shown that an arbitrary element x є V can be expressed uniquely in the form 


Х = Хуур + + X,Y, (x; € F). 


In physical and engineering applications, it is useful to distinguish between physical 
quantities which are associated with a direction such as velocity, acceleration, force 
and magnetic field, and quantities with no associated direction such as mass, density 
and temperature. The latter can be represented by real numbers, while geometric 
vectors may be used to handle the former. 


= 
Geometric vectors are equivalence classes of parallel arrows of equal length. If AB is 
an arrow, we write AB for the geometric vector to which it belongs. The sum of two 


geometric vectors is obtained by choosing suitable arrows which can be added accord- 
ing to the triangle law of addition: 


OA + AB = OB. 
> => — 


о 


If 4 is a real number then we define the product 2AB to be the geometric vector 
whose arrows are parallel to those of AB but with lengths multiplied by 2. 
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With respect to these definitions of addition and scalar multiplication the set G? of 
geometric vectors in space forms a real vector space. We often consider the set G? of 
geometric vectors in the plane, which is a subspace of this vector space. The null 
vector 0 is the geometric vector whose arrows are all of zero length. 


Let Oxyz be a frame of rectangular Cartesian coordinate axes and let i, j, k be 
geometric vectors whose arrows are of unit length and parallel to the coordinate axes 
Ox, Oy, Oz respectively. Then every geometric vector in the plane can be written 
uniquely in the form 


а= аі + aj, 
and every geometric vector in space can be written uniquely in the form 
a= aji + a,j + ask. 


Thus fi, j} forms a basis for the space С? of geometric vectors in the plane, and (i, j, k} 
forms a basis for G?. a,, a5, ау are, respectively, the x-, y-, z-coordinates of a. 


Another important vector space is R?, the set of ordered pairs of real numbers, with 
addition and scalar multiplication defined by 


(a1, 2) + (b, b2) = (а, + by, a2 + ba), 
and 
m(a,, a5) = (ma, , maz). 


R? isa real vector space with respect to these operations and we frequently identify R? 
with G? by the correspondence 


a,i + аз) e» (a1, 4), 


which is an isomorphism. Similarly we may identify R? with G?. Under these identifica- 
tions, OA (where O is the origin of coordinates) corresponds to the element of R? 


or R? which gives the coordinates of A in the coordinate system of Oxyz. 


A 


Cler со. з) 


A(a,. 2, 43) 


ау 


Bibi. ba. bs.) 


x 


— — 
If OA and BC belong to the same geometric vector then 


а, = сү — by, a, = с; — bz, a3 = сз — bs, 
where 

OA = aj a,j + ak, 

OB = bj + bj + bk, 

OF = сі + сј + c3k. 


OA or (a,, 25, a3) is the position vector of A. 
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An important property of a geometric vector is its length or magnitude. To define this 
more generally, we introduce the inner product. 


A mapping -: V x V — R, where V is а real vector space, is called a (real) inner 
product if for all a, b, ce Vand A, ue R: 


1 а:а= 0а = 0; 

2 a:a20; 

3 a:b-b-a; 

4 a-(4b + pe) = да: + pase. 

We may specify an inner product on G? or G? by 
ОА. ОВ = OA. OB cos 0, 


— —— 
where OA, OB denotes the length of thearrow OA, OB respectively and 0 is the angle 
— — 
between OA and OB. The geometric vectors i, j, k are orthogonal, i.e. 


је р к=к = 0. 


(In fact, any two geometric vectors which are represented by perpendicular arrows аге 
orthogonal.) Thus 


OA ОВ = a,b, + azb, + а,Ь; 

in terms of coordinates. Correspondingly, we define, for vectors in R?, 
(a), a2, аз) * (by, bz, b3) = a,b, + aiba  asbs. 

This satisfies the axioms for an inner product on R?. 

The norm of a vector is defined by 
lal = (а ғ a}. 


Thus the norm ofa geometric vector is just the length of an arrow representing it, and 
the norm of an element of R? is given by 


lai, а, az) = (af + a} + a3). 


We adopt the convention that a denotes |||, i.e. if a vector is denoted by a boldface 
symbol then the corresponding italic symbol denotes its norm. Note that a* a = a’. 


An important theorem in real inner product spaces is Schwarz's Inequality: 
(xy)? € (x: x)(y- y). 

In G? (but not in С?) we define the vector product x :G? x G? — G? by 
OA x OB = OA.OBsin 0n. 


— — 
where 0 is the angle (0 < 0 < x) between QA and OB, and nis the unit vector (i.e. vector 


whose norm is 1) orthogonal to ОА and OB, such that OA, OB and n form a right- 
handed triad. 
The vector product has the following properties: 
1 ixj=kjxk=ikxi=jixi=jxj=kxk=0; 
2 axb--bxa; 
3 ax(b+c)=axbt+axc; ; 
4 axbc-(ajbs — а,Ь») + (azb; — ayb3)j + (aiba — ajbj)k; 
5 (a x b) xc #a x (bxc), in general. 
Property 4 may be obtained by expanding 
a x b = (aii + a5j + ask) x (bib + bj + БК), 


using Properties 1, 2 and 3 which are simple consequences of the definition. 
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If ABCD ‘is a parallelogram, i.e. 
A= DG and AD = BG 
then its area is given by 
[АВ х BC]. 
— — 


The scalar obtained by evaluating a : (b x c) is called the triple scalar product of 
2, b and c. For convenience, we often write a -b x c; no confusion should arise since 
(a:b) x c would be meaningless. We note that 


а:Ьхе= |а b, с, 
аз b} C3}, б 
and, by the property of determinants regarding interchange of columns, we obtain 
a‘bxc=b'cxa=c:axb. 
The triple scalar product determines a function from G? x G? х G? — R. 


We may show that the volume of a parallelepiped (i.e. a three-dimensional solid whose 
faces are parallelograms) is given by the modulus of the triple scalar product 


la-b x el, 


where a, b and c are geometric vectors represented by three arrows which generate the 
parallelepiped. 
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PARTIAL DIFFERENTIAL EQUATIONS OF APPLIED MATHEMATICS 


1 W The Wave Equation 
2 W Classification and Characteristics 


3 W Elliptic and Parabolic Equations 


4 NO TEXT 

5 S Finite-Difference Methods I: Initial Value Problems 
6 W Fourier Series 

7 N Motion of Overhead Electric Train Wires 

8 S Finite-Difference Methods II: Stability D 
9 W Green's Functions 1: Ordinary Differential Equations 
10 W Green's Functions II: Partial Differential Equations 


11 S Finite-Difference Methods III: Boundary Value Problems 
12 NO TEXT 

13 W Sturm-Liouville Theory 

14 W Bessel Functions 

15 N Finite-Difference Methods IV 

16 N Blood Flow in Arteries. 


The letter after the unit number indicates the relevant set book; N indicates a unit 
not based on either book. 
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